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strain  rates  of  5  x  10’  s“^  .v^The  next  section  describes  the  results  of  a 
detailed  study  of  the  decomposition  of  PETN  using  time-of-flight  mass  spec¬ 
trometry^  The  decomposition  was  induced  by  heating  slowly,  by  fracturing 
single  crystals  and  by  laser  irradiation. ^  Different  reaction  pathways  were 
found  in  each  case  and  these  are  discussed. ^Finally,  an  experimental  study 
is  described  in  which  the  deformation  stress  of  millimeter-size  single 
crystals  of  explosive  materials  has  been  determined  using  a  novel  experi¬ 
mental  technique^  The  method  involved  the  loading  of  the  test  crystal 
against  a  transparent  glass  plate  with  in  situ  measurements  of  the  real 
area  ot  c on tacy  under  load  using  an  optical  microscope.  For  all  explosive 
crystals  the  deformation  stress  was  found  to  increase  with  load  reaching  a 
plateau  at  higher  loads.  The  maximum  deformation  stress  of  a  crystal  was 
found  to  be  considerably  lower  than  its  Vickers  diamond  hardness  value. 
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DEFORMATION,  FRACTURE  AND  EXPLOSIVE  PROPERTIES 
OF  REACTIVE  MATERIALS 


1.  INTRODUCTION 

In  this  annual  report,  we  describe  four  areas  of  work  completed  during 
the  grant  period.  This  research  has  recently  been  prepared  for  publication 
and  the  full  texts  of  the  papers  are  given. 

1.1  Drop^welght  impact  studies 

These  have  been  made  using  the  facilities  developed  in  early  work  and 
described  by  Heavens  and  Field  (1974).  Figure  1,  of  the  first  paper, 
illustrates  our  transparent  anvil  drop-weight  apparatus  which  allows  us  to 
photograph  events  at  microsecond  framing  rates.  Figure  2  is  an 
instrumented  drop-weight  machine  which  allows  us  to  record  pressure-time 
curves . 

In  the  grant  period,  we  have  made  a  detailed  study  of  hexanitrostllbene 
(HNS)  and  have  compared  its  behaviour  with  a  range  of  other  explosives. 
The  photographic  sequences  show  that  upon  impact  the  layer  of  explosive 
undergoes  bulk  plastic  flow  at  speeds  of  over  250  m  s_1 .  Initiation  occurs 
after  this  catastrophic  failure.  From  a  comparison  of  melting  points, 
uniaxial  yield  stress  values,  flow  properties  and  the  ignition 
characteristics  of  several  organic  secondary  explosives,  it  appears  that 
explosives  with  high  strength  values  are  more  likely  to  exhibit  both  a 
rapid  mechanical  failure  and  resultant  ignition  during  impact.  The  results 
emphasize  the  important  relation  between  mechanical  properties  and  ignition 
characteristics. 

1 .2  High  strain  rate  properties  of  explosives 


This  second  paper  was  presented  at  the  Second  International  Conference 


on  "The  Mechanical  Properties  of  Materials  at  High  Rates  of  Strain",  held 
at  Oxford,  U.K.  in  April  198*4.  The  paper  describes  our  experimental 
facilities  which  include:  (a)  drop-weight  machines;  (b)  a  technique  which 
uses  the  so-called  "Brazilian"  test  geometry  and  laser  speckle  photograph/ 
for  measurements  of  the  tensile  strengths  and  rupture  strains  of  PBX 
explosives;  and  (c)  a  miniaturised  and  modified  Hopkinson  bar  apparatus 
which  allows  us  to  record  the  stress  strain  curves  of  explosives  at  strain 
rates  up  to  10^  s-1 .  Results  are  given  in  this  paper  for  various  A.W.R.E. 
Aldermaston  PBX  compositions.  Recently  we  have  started  research  on  U.S. 
compositions  supplied  by  Los  Alamos. 

1 .3  Thermal ,  fracture  and  laser-induced  decomposition  of  PETN 

This  project  which  involved  a  detailed  mass  spectrometry  study  was 
started  a  few  years  ago,  but  has  recently  been  completed  and  prepared  for 
publication.  It  describes  the  different  reaction  pathways  following  slow 
heating,  fracture  and  laser-irradiation  of  PETN. 

1 . *4  The  deformation  stress  of  highly  brittle  explosive  crystals  from  real 

contact  area  measurements 

An  experimental  study  is  described  in  which  the  deformation  stress  of 
mm.  size  highly  brittle  as-grown  single  crystals  of  explosive  materials  has 
been  determined  using  a  novel  experimental  technique.  The  technique 
Involved  the  loading  of  the  test  crystal  against  a  transparent  glass  plate 
with  in  situ  measurements  of  the  real  area  of  contact  under  load  using  an 
optical  microscope.  For  all  explosive  crystals  the  deformation  stress  was 
found  to  increase  with  load  reaching  a  plateau  at  higher  loads.  The 
maximum  deformation  stress  of  a  crystal  was  found  to  be  considerably  lower 
than  its  Vickers  diamond  hardness  value.  A  discussion  on  the  determination 


of  the  uni-axial  yield  stress  of  small  particles  is  given. 
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2.  IMPACT  INITIATION  OF  HEXANITROSTILBENE 
V.  Krishna  Mohan  and  J.E.  Field 


ABSTRACT 

A  study  has  been  made  of  the  Initiation  of  hexanitrostilbene  (HNS) 
under  impact  using  high-speed  photography  at  microsecond  framing  rates. 
The  photographic  sequences  show  that  upon  impact  the  layer  of  explosive 
undergoes  bulk  plastic  flow  at  speeds  of  over  250  m  s“1.  Initiation  occurs 
after  this  catastropic  failure  of  the  explosive  due  to  the  fact  that  the 
adiabatic  heating  produced  by  plastic  deformation  is  extremely  localized  in 
nature.  From  a  comparison  of  the  melting  points,  uniaxial  yield  stress 
values,  flow  properties  and  the  ignition  characteristics  for  several 
organic  secondary  explosives,  it  appears  that  explosives  with  high  strength 
values  are  more  likely  to  exhibit  both  a  rapid  mechanical  failure  and 
resultant  ignition  during  impact.  Addition  of  polycarbonate  (PC)  has  been 
found  to  sensitize  the  explosive.  Strain-gauge  measurements  of  impact  of 
pure  HNS  and  HNS/PC  samples  have  shown  that  a  sudden  pressure  drop  occurs 
before  ignition  due  to  rapid  failure  of  the  test  material,  the  latter  being 
in  agreement  with  the  photographic  observations. 


4. 


.  INTRODUCTION 

Hexanitrostilbene  (HNS)  is  a  secondary  explosive  which  is  used  in  many 
applications  because  of  its  relative  stability  at  elevated  temperatures. 
Of  major  concern  to  explosive  formulators  is  the  possibility  of  accidental 
initiation  during  manufacture,  transport  or  usage.  One  of  the  standard 
methods  adopted  to  assess  the  hazards  associated  with  an  explosive  is  the 
fall-hammer  test  in  which  a  sample  (typically  20  -  50  mg)  is  subjected  to 
impact  with  a  falling  weight.  Repeated  series  of  these  tests  give  some  idea 
about  the  impact  energy  necessary  to  ignite  the  sample.  The  results  of 
these  tests  are  usually  expressed  in  terms  of  H50*  where  H50  refers  to  the 
height  from  which  the  weight  would  ignite  50%  of  a  series  of  samples. 
Although  the  test  has  been  used  for  many  decades,  it  is  only  recently  that  a 
systematic  attempt  has  been  made  to  photograph  the  sample  during  its 
deformation  [1*3] .  A  combination  of  high-speed  photographic  and 

pressure-measuring  techniques  helped  establish  the  physical  processes 
occurring  during  Ignition  and  propagation. 

The  high-speed  photographic  studies  pioneered  by  Heavens  and  Field  [1] 
have  been  useful  in  identifying  the  various  mechanisms  responsible  for  the 
formation  of  initiation  sites.  In  general  it  is  thought  that  the  initiation 
of  explosion  by  mechanical  shock  is  thermal  in  origin.  On  the  basis  of 
localized  thermal  energy  or  'hot-spots'  as  the  source  of  explosion,  four 
possible  mechanisms  for  ignition  by  impact  can  be  envisaged: 

(i)  adiabatic  compression  of  trapped  gas  spaces  [4,5] 

(li)  viscous  heating  of  material  rapidly  extruded  between  the  impacting 
surfaces  or  grains  [6,7] 

(ili)  friction  between  the  impacting  surfaces,  the  explosive  crystals 
and/or  grit  particles  in  the  explosive  layer  [8,9] 

(lv)  localized  adiabatic  shear  of  the  material  during  mechanical  failure 

[3,  10  -  12] 

The  references  cited  above  for  the  various  initiation  mechanisms  are  not 


exhaustive;  a  complete  list  is  given  in  reference  3*  Field  et  al  [3]  have 
shown  that  no  one  mechanism  is  the  dominant  means  of  ignition  and  that  only 
very  slight  changes  in  the  conditions  of  an  experiment  can  lead  to 
•hot-spots'  being  formed  by  different  processes. 

The  high-speed  camera  work  of  Heavens  and  Field  [1]  revealed  that 
various  phenomena  may  be  involved  when  a  thin  layer  of  material  is  impacted, 
including  plastic  flow,  jetting,  shear  banding  and  even  melting.  Ignition 
occurred  at  a  few  'hot-spots*  following  severe  deformation  of  the  sample. 
Strain-gauge  measurements  showed  that  if  a  sample  failed  by  plastic  flow,  a 
sharp  pressure  drop  occurred  of  a  few  hundred  MPa.  It  was  after  the 
pressure  drop  that  ignition  took  place.  Of  the  several  explosives  examined 
by  these  authors,  only  pentaerythritol  tetranitrate  (PETN), 
cyclotrimethylene  trinitraraine  (RDX)  and  cyclotetramethylene  tetranitramine 
(HMX )  exhibited  such  large  pressure  drops  and  ignited  after  the  pressure 
drop  while  other  explosives  such  as  tetryl  and  picric  acid  gave  rise  to  a 
series  of  pressure  fluctuations  to  which  ignition  was  not  linked.  Recent 
photographic  studies  on  the  impact  ignition  of  tetryl  indicated  that  it  has 
comparatively  little  tendency  to  flow  under  impact;  the  explosion  was 
quite  localized  at  the  various  ignition  sites  close  to  the  periphery  of  the 
sample  and  did  not  propagate  into  the  bulk  of  the  sample  [13]. 

It  is  clear  from  the  above  discussion  that  there  are  differences  in  the 
mode  of  ignition  for  different  explosives  and  It  is,  therefore,  worthwhile 
to  widen  the  data  base  by  examination  of  further  explosive  materials.  With 
this  objective  in  view,  the  impact  behaviour  of  HNS  is  investigated  in  the 
present  work  using  both  photographic  and  pressure-measuring  techniques. 

2.  EXPERIMENTAL 

High-speed  photographic  observation  of  the  behaviour  of  a  thin  layer  of 
explosive  was  carried  out  using  toughened  glass  anvils,  allowing  the  impact 
process  to  be  viewed  in  transmitted  light.  An  A.W.R.E.  designed  C  4 
rotating-mirror  camera  at  framing  Intervals  of  ca.  7.0  us  was  used  to 


photograph  the  event.  Figure  1  shows  the  experimental  arrangement  of  the 
impact  set-up;  details  of  the  apparatus  are  given  in  -’eference  1 .  In  this 
work,  all  experiments  were  made  on  25  mg  of  HNS  in  crystalline  form 
(particle  size  50  -  150  ym),  which  was  impacted  by  a  5.5  kg  weight  falling 
from  a  height  of  1  m.  A  few  trials  were  carried  out  by  adding  polycarbonate 
(PC)  to  the  explosive. 

Experiments  were  also  performed  on  an  instrumented  drop-weight  machine 
similar  to  that  described  by  Heavens  and  Field  [1]  (see  Figure  2).  The 
apparatus  was  connected  to  a  self  contained  data  recording  system.  This 
comprised  a  Data  Laboratories  D1922  transient  recorder,  a  Tektronix 
double-beam  oscilloscope,  and  a  paper  tape  punch.  With  this  system, 
analogue  signals  from  the  drop-weight  apparatus  are  digitised  in  real  time 
and  retained  in  memory  30  that  the  signals  may  be  subsequently  photographed 
on  the  oscilloscope  and  transferred  to  the  paper  tape  for  data  analysis. 
The  data  was  analysed  on  a  HP9825  calculator  by  using  a  program  based  on  the 
method  or  analysis  outlined  in  references  1  and  10. 

3.  RESULTS  AND  DISCUSSION 

Photographic  observation  of  the  Impact. 

A  sequence  illustrating  the  behaviour  of  HNS  under  impact  is  given  in 
Fig.  3-  Frame  (a)  corresponds  to  the  initial  instant  of  impact  while  frames 
(b)  and  (c)  respresent  the  sample  after  315  and  1*34  ys  respectively.  During 
the  early  stages  the  thin  explosive  layer  is  compacted  into  the  form  of  a 
pellet.  By  frame  (d)  the  layer  starts  to  deform  plastically  and  from  this 
stage  onwards  the  flow  of  material  is  extremely  rapid,  attaining  speeds  of 
over  250  m  s-1.  It  appears  that  the  explosive,  initially  opaque,  becomes 
translucent  as  it  undergoes  severe  plastic  deformation  presumably  due  to 
fusion  of  particles  under  impact  pressure.  Ignition  occurs  after  ca.  450  y 
at  the  periphery  of  the  expanding  layer  and  subsequently  propagates  int< 
the  bulk  of  the  sample.  The  later  frames,  (f)  to  (i)  show  that  the  centra] 
portion  which  is  translucent  in  frame  (e)  gradually  becomes  opaque  as  the 


3.  HIGH  STRAIN  HATH  I'llUTLllUES  OF  EXPLOSIVES 


J.E.  Field,  G.M.  Swallowe,  P.II.  Pope  and  S.J.P.  Palmer 
Cavendish  Laboratory,  Madingley  Hoad,  Cambridge  CB3  OilE 


Abstract..  Experiments  are  described  on  the  deformation  and  iGnition 
behaviour  of  a  range  of  explosives,  including  some  polymer  bonded 
explosives  (PCX's).  An  instrumented  drop-weight  apparatus  wan  used  to 
obtain  stress-strain  curves  for  strain  rates  in  tho  range  I02  to  1 0 3  s“1. 
A  novel  drop-weight  apparatus  with  transparent  anvils  allowed  high-speed 
photographic  recording  of  the  impact  and  ignition  processes.  Laser 
speckle,  used  in  conjunction  wi  th  a  specimen  loaded  In  the  Brazilian 
test  geometry  gave  tensile  strength  and  strain  to  failure  data  of  PBX'r. 
Finally,  a  r.odlfied  llopkinson  bar  apparatus  gavo  data  at  strain  rates  up 
to  5  x  101*  s-'. 

1.  Introduction 


It  is  generally  accepted  that  explosive  initiation  is  thermal  in  origin.  Jn 
processes  involving  impact  or  shock,  the  mechanical  energy  is  envisaged  ns 
being  converted  into  heat  in  localized  regions  called  "hot  spots"  (Bowden 
and  Yoffe  1952, 1 9tiC)  the  required  "hot  spot"  radius  for  ignition  is  a 
function  of  the  explosive,  the  temperature  arid  its  duration  (Rideal  and 
Robertson,  1 9 *>8 ) .  For  pentaerythritol  tetranitrrtc  (PETR),  the  calculated 
"hot  spot"  temperature  for  a  time  of  1  ps  and  radius  1  pin  is  ca.  700  K. 
Experimental  work  by  Bowden  and  Gurton  (19*i9a,b)  largely  confirmed  these 
figures  and  showed  that  "hot  spot"  radii  are  typically  0.1  to  10  pm  for 
Initiation  by  friction  or  impact.  In  a  recent  paper  (Field  ct  al ,  198?.)  wo 
have  discussed  the  ignition  mechanisms  of  explosives  during  mechanical 
deformation.  Evidence  was  given  for  ignition  by  adiabatic  heating  of 
trapped  gas,  adiabatic  shear  of  the  explosive,  friction,  viscous  flow, 
fracture  or  shear  of  added  particles  and  tri bolumineoeent  discharge. 

A  feature  of  our  work  on  the  relation  between  mechanical  properties  anil 
explosive  behaviour  has  been  the  use  of  high-speed  photography.  The 
advantage  of  the  piiotographi c  evidence  is  that  it  has  allowed  us  to  identify 
the  conditions  under  which  the  various  ignition  mechanisms  oporate  during 
deformation  of  ari  '  (plosive. 

In  many  modern  applications,  explosive  crystals  are  bonded  into  a  polymer 
matrix  (so-called  PbX  explosives).  It  is  important  with  sucli  systems  to 
understand  the  factors  which  affect  their  mechanical  and  therm:*! 
properties,  their  sensitiveness  (response  to  a  prescribed  stimulus)  and 
their  cxplosi  vfiwr.r,  ((/plosive  response).  There  is  usually  a  conflict 
between  these  various  requirements  arid  an  optimum  has  to  he  soucht.  Tin 
present  paper  describes  results  on  pure  explosives,  explosives  with 
additives  and  some  f’HX  explosives  based  on  eye!  of  el .  rain'd,  hy  1  <nc 
tetrani  Irani  no  (II, MX  ) , 
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Figure  1  Experimental  arrangement  at  instant  of  impact.  W,  drop 

weight;  G,  glass  blocks;  M,  mirror;  S,  sample.  The  upper 
glass  block  G  is  attached  to  the  drop  weight. 

Figure  2  Experimental  arrangement  for  obtaining  pressure-time  curves. 

W,  drop-weight;  H,  R12,  R2,  R3  hard  steel  rollers; 

C,  cylindrical  guiding  sleeve;  S,  sample;  G,  strain  gauge; 
A,  cast  steel  anvil. 

Figure  3  Impact  on  a  layer  of  granular  HNS.  The  sample  becomes 

translucent  [frame  (d ) ]  and  starts  to  undergo  plastic 
flow  [frame  (e)].  Ignition  occurs  in  frame  (f)  at  the 
periphery  (see  arrows).  From  frame  (g)  onwards  fast  reaction 
is  taking  place  in  the  bright  areas  indicated  by  the  arrows 
(3ee  frame  h). Frame  tim<'/ys  :  (a)  0;  (b)  315; 

(c)  434;  (d)  441;  (e)  448;  (f)  455;  (g)  490;  (h)  518 
(i)  560. 

Figure  4  Photograph  of  the  toughened  glass  anvils  after  the  impact 

test  with  (a)  pure  explosive  (unreacted  material  indicated 
by  an  arrow)  and  (b)  explosive/PC  sample. 

Figure  5  Radius  of  a  layer  of  HNS  during  Impact  showing  evidence 

of  plastic  flow. 

Figure  6  Stress  time  curves  for  impact  of  (a)  HNS  and  (b)  HNS/PC. 

Mass  of  hammer  5  kg;  drop  height  1  m.  The  initially 
lower  trace  records  the  instant  of  ignition  by  monitoring 
the  change  in  electrical  resistance  of  the  sample.  The 
trace,  labelled  S,  is  the  stress-time  curve. 

Figure  7  Pressure-time  curves  for  impact  of  (a)  HNS  and  (b)  HNS/PC. 


TABLE  2 


Strength  properties  and  impact  behaviour 
of  various  materials. 


Compound 

Melting 

point/K 

Yield 

stress/MPa 

H50/m  Observations  in 
Rapid  flow 

impact  apparatus 

Ignition 

TNT 

35  4a 

34b 

1.48° 

Nod 

No 

Picric  acid 

395 

52 

0.73 

No 

No 

Tetryl 

40^4 

53 

0.37 

No 

Yese 

PETN 

414 

60 

0.13 

Yes 

Yes 

RDX 

474 

82 

0.28 

Yes 

Yesf 

Ammonium 

perchlorate 

>493 

— 

— 

Yes 

Yes 

HMX 

553 

125 

0.32 

Yes 

Yes 

HNS 

588 

-  140 

0.54 

Yes 

Yes 

Potassium 

nitrate 

>660 

Yes 

a.  Values  taken  from  Lawrence  Livermore  National  Laboratory 
Explosives  Handbook,  UCRL-52997  (1981). 

b.  All  the  yield  stress  values  from  ref.  [10], 

except  for  HNS  which  has  been  estimated  using  equation  (3). 

c.  Drop-weight  2.5  kg;  sand  paper  surface.  Data  from  reference 
given  above  in  "a". 

d.  The  experimental  observations  for  all  compounds  excluding 
tetryl  and  HNS  are  from  ref.  [1]. 

e.  Reaction  fails  to  propagate. 

f.  Ignition  with  large  masses  (100  mg)  only. 
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TABLE  1 


Infra-red  absorption  frequencies  and 
corresponding  assignments  for  Hexanitrostilbene 


Frequency/cm-1 

Assignments 

725 

(3) 

Aryl-CH-CH-Aryl  out  of  plane  deformation. 

740 

(s) 

N0g  out  of  plane  deformation;  inconsistent 

band,  limited  value. 

810 

(a) 

Aromatic  (Ar)C-H  out  of  plane  deformation. 

825 

(a) 

C-N  Vibration;  limited  value. 

920 

(a) 

Ar  C-H  out  of  plane  deformation. 

940 

(a) 

955 

(a) 

Aryl-CH-CH-Aryl  out  of  plane  deformation. 

1080 

(m) 

Ar  C-H  in  plane  deformation. 

1340 

(a) 

Symmetric  NO2  stretching. 

1460 

(b) 

Ar  C-C  in  plane  deformation. 

1535 

(a) 

Asymmetric  NOg  stretching. 

1600 

(a) 

Ar  C-C  in  plane  deformation. 

1620 

(m) 

Ar  C-C  stretching. 

3100 

(a) 

-  C-H  streching. 

Letters  in  parenthesis  refer  to  the  intensity  of  the  band:  s,  sharp; 
m,  medium;  b,  broad. 
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explosives  which  possess  high  strength  are  more  susceptible  to  catastropic 
failure  and  ignition  during  impact.  This  helps  explain  the  behaviour  of 
different  materials  in  the  drop-weight  test  and  has  important  practical 
consequences  for  explosive  formulations. 
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cubic  metals  and  semiconductors  shows  a  qualitative  correlation  between 
these  two  quantities  [17].  These  facts  suggest  a  link  between  the  elastic 
moduli  of  a  material,  whose  magnitude  is  to  a  large  extent  determined  by  the 
strength  of  the  inter-atomic  (or  molecular)  bonding  forces,  and  the  melting 
temperature  which  is  another,  albeit  indirect,  measure  of  such  forces. 
Equation  (3)  is  not  valid  in  the  case  of  secondary  explosives  because  the 
units  of  crystal  structures  in  such  compounds  are  comprised  of  complex 
molecules,  which,  in  turn,  exhibit  a  vibrational  complexity.  However,  for 
such  secondaries,  the  uniaxial  yield  stress,  which  is  another  indicator  of 
the  mechanical  strength  of  a  solid  is  found  to  increase  with  the  melting 
temperature  of  the  explosive  [10] (see  Table  2).  A  similar  trend  is  noticed 
in  the  hardness  data  for  some  of  these  explosives.  For  example,  the  Vickers 
hardness  values,  in  units  of  kg/mm^,  for  PETN,  on  it3  cleavage  plane  {110}, 
RDX,  on  the  (210)  growth  surface,  and  B-HMX,  on  its  cleavage  plane  {Oil}, 
are  17.9  (ref.  18),  35.2  (average  of  values  given  in  ref.  19)  and  HO. 3  (ref. 
20),  respectively,  while  the  corresponding  melting  points  are  HlH,  H7H  and 
533  K-  The  hardness  value  for  HNS  is  not  known  but  it  is  reasonable  to 
expect  that  it  will  be  high.  Young's  modulus,  uniaxial  yield  stress  and 
Vickers  hardness  are  related  to  each  other  by  different  expressions,  see 
for  example  ref.  [20].  Table  2  also  gives  the  yield  stress  values  and  the 
impact  characteristics  such  as  flow  properties  and  initiation  of  fast 
reaction  for  a  number  of  explosives  and  also  for  potassium  nitrate  (which  is 
a  non-explosive  compound).  These  data  suggest  that  the  strength  of  a 
compound  in  the  drop-weight  test  is  a  key  factor  in  determining  the  impact 
behaviour.  If  the  strength  is  low,  then  deformation  starts  at  a  low  stress 
and  flow  proceeds  at  moderate  velocities.  On  the  other  if  it  is  high,  the 
onset  of  deformation  is  delayed  but  then  develops  catastrophically. 
Examples  of  flow  at  many  hundreds  of  m  s-1  have  been  recorded  in  our  work. 
Two  important  factors  in  the  generation  of  'hot-spots'  are  the  rate  of 
deformation  and  the  resistance  to  this  deformation.  Thus  it  appears  that 
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prerequisite  for  ignition  is  in  agreement  with  those  of  Kholevo  [15], 
Afanasev  and  Bobolev  [10]  and  Heavens  and  Field  [1],  As  mentioned  earlier 
the  failure  of  the  sample  is  rapid  in  nature. 

Strength  and  sensitivity  relation 

A  comparison  of  the  flow  properties  for  several  explosives  based  on  the 
results  of  Heavens  and  Field  [1]  and  our  work  on  HNS  and  tetryl  [13] 
suggests  that  among  organic  secondary  explosives,  those  with  a  high  melting 
point  normally  exhibit  a  greater  tendency  to  flow  rapidly  under  impact. 
Thus  trinitrotoluene  (TNT),  picric  acid,  and  tetryl  with  melting  points  in 
the  range  353-i*0J4  K  show  only  minor  deformation  upon  impact  while 
explosives  with  high  melting  points  such  as  RDX  and  HMX  undergo  bulk  plastic 
flow  in  the  glass  impact  apparatus.  HNS  is  a  high  melting  compound  (melting 
temperature  588  K)  therefore,  the  extensive  flow  seen  in  impact  tests 
appears  in  line  with  previous  results.  The  borderline  for  rapid  flow  in 
this  equipment  appears  to  lie  between  tetryl  and  PETN.  The  only  other 
explosive  studied  which  melts  throughout  the  bulk  is  RDX,  but  then  only 
between  steel  anvils  [1,3]. 

For  ionic  solids  with  a  simple  cubic  structure  a  semi-empirical 
relation,  based  on  the  Linderaann  theory  of  vibrational  instability  of 
crystal  lattices,  can  be  derived  between  the  melting  temperature  T  and  the 
Young's  modulus  E 

T  -  Hr- 62  (3) 

where  p  is  the  density  of  the  solid  in  kg/m3,  m  is  the  molar  mass  in  kg  and 
R  is  the  gas  constant  in  J/mol/K  [16].  The  vibration-amplitude  x0  when 
melting  occurs  is  related  to  the  atomic  (or  molecular)  spacing,  a,  by 
xc  -  6a  with  0  <  6  <  1.  Although  equation  (3)  is  based  on  a  crude  model,  it 
has  yielded  good  agreement  between  the  calculated  and  the  observed  melting 
points  for  several  materials.  Also,  a  plot  of  E  versus  T  for  a  number  of 


impact  apparatus,  the  reaction  always  proceeds  to  completion  in  this  test. 
This  is  because  the  modulus  and  hardness  of  the  steel  used  are  greater  than 
that  of  toughened  glass.  With  HNS/PC  two  small  incremental  drops  were 
recorded  soon  after  impact  and  are  attributed  to  the  rapid  deformation  of 
the  polymer. 

The  impact  presure,  P,  is  given  by  the  relation 

P  -  aU  (1) 

where  the  constant  a  is  the  gauge  calibration  factor  and  U  is  the  voltage 
signal  from  the  gauge.  The  gauge  was  calibrated  statically  using  an  Instron 
machine  and  the  calibration  factor  was  found  to  be  0.043  GPa/mV.  The 
pressure-time  curves  obtained  using  equation  (1)  are  presented  in  Fig.  7. 
At  first  the  sample  is  compressed  and  compacted.  Eventually  the  strength  of 
the  layer  is  exceeded  and  the  material  deforms  plastically.  In  the 
drop-weight  test  the  average  pressure  for  yield,  Puit*  is  higher  than  the 
uniaxial  yield  stress  of  the  material  oy.  This  is  due  to  the  constraining 
effect  of  the  anvils.  Application  of  the  Mises  yield  criterion  to  a  thin 
layer  of  thickness  h  and  diameter  d  between  two  rigid  cylindrical  punihes 
gives  the  following  hyperbolic  expression  [14], 

Pult  -  0y[1  +  (d/3/3  h)]  (2) 

Note  that,  the  smaller  the  thickness,  or  greater  the  diameter,  the  higher 
the  value  of  Puit*  From  Fig.  7  it  can  be  seen  that  Pu^t  value  for  HNS  under 
the  test  conditions  is  760  MPa.  For  HNS/PC  two  pressure  drops  are  recorded; 
the  lower  values  for  the  explosive/polymer  sample  suggesting  that  they 
correspond  to  the  failure  of  PC.  Using  equation  (2),  an  approximate  value 
for  the  yield  stress  of  HNS  can  be  estimated.  oy  is  thus  found  to  be  140 
MPa. 

For  HNS  both  the  impact  photographs  and  the  stress-time  curves  show 
that  there  is  a  close  relation  between  the  mechanical  failure  of  the  layer 
and  the  instant  of  ignition;  the  ignition  always  occurring  either  during  or 
after  the  failure.  Our  observation  that  rapid  deformation  is  a  necessary 
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under  impact  is  'partial1.  Swallowe  and  Field  [2]  have  examined  in  detail 
the  effect  of  a  number  of  polymers  on  the  H50  values  of  the  secondary 
explosives,  PETN  and  HMX.  These  aut'  'Ts  found  that  some  relatively  soft, 
low-melting  polymers  such  as  PC  and  polysulphone  (PSu)  can  also  sensitize 
explosives.  These  polymers  fail  catastrophically  either  by  fracture  or 
localized  adiabatic  shear  and  have  a  low  specific  heat,  latent  heat  and 
thermal  conductivity.  The  temperature  of  the  'hot'spots'  in  such  polymers, 
produced  during  rapid  deformation,  can  greatly  exceed  the  polymer's 
softening  point.  Neither  the  size  of  the  grit  nor  the  percentage  content  in 
the  explosive/grit  mixture  were  seen  to  influence  the  extent  of 
sensitization.  For  HMX,  addition  of  PC  reduced  the  H50  values  from  29.0  ± 
2.0  cm  (for  the  pure  compound)  to  16. 0  t  1.5  cm  [2].  The  sensitization  is 
due  to  'hot  spots'  generated  by  non-homogeneous  energy  dissipation  either 
in  adiabatic  shear  bands  or  at  crack  tips  in  the  polymer.  In  the  current 
work  the  effect  of  PC  on  the  impact  initiation  of  HNS  was  investigated.  In 
three  trials  with  HNS/PC  samples  between  glass  anvils,  explosion  always 
propagated  to  completion  with  no  explosive  remaining  on  the  anvils.  The 
bottom  anvil  from  one  such  trial  is  shown  in  Fig.  4(b).  In  the  steel  impact 
apparatus,  the  H50  value  of  the  pure  explosive  decreased  from  0.68  ±  0.02  m 
to  0.40  ±  0.03  m  in  the  presence  of  the  polymer. 

Impact  behaviour  in  the  drop-weight  test 

The  stress-time  curves  for  pure  HNS  and  HNS/PC  were  obtained  using  an 
instrumented  drop-weight  machine  with  steel  anvils  and  are  presented  in 
Fig.  6(a)  and  (b)  respectively.  The  trace  (marked  S)  represents  the 
stress-time  relation  recorded  by  the  strain  gauge  during  impact  and  the 
lower  trace  detects  the  Instant  of  ignition  by  recording  a  jump  in  the 
electrical  conductivity  of  the  material  between  the  anvils.  For  the  pure 
explosive,  the  prominent  feature  of  Fig.  6  (a)  is  a  sharp  break  in  the 
stress-time  curve  associated  with  ignition  following  mechanical  failure  of 
the  sample.  Contrary  to  the  'partial'  native  of  reaction  in  the  glass 
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explosion  propagates  further.  The  transition  from  a  translucent  to  an 
opaque  layer  is  thought  to  be  due  to  the  fact  that  the  gaseous  explosion 
products  separate  the  impacting  anvils  slightly  and  reduce  the  impact 
pressure.  The  photographs  also  show  that  the  fast  reaction  does  not 
propagate  into  the  centre  of  the  sample.  In  fact,  examination  of  the  bottom 
anvil  after  the  test  indicates  that  a  certain  amount  of  explosive  remains 
behind.  The  anvil  is  shown  in  Fig.  Ma).  In  order  to  ascertain  that  the 
residue  is  unreacted  explosive  rather  than  partially  decomposed  material, 
its  infra-red  spectrum  was  taken  in  Nujol  and  compared  with  that  of  pure 
HNS.  The  two  spectra  were  found  to  be  indentical  thus  confirming  that  the 
material  on  the  anvil  was  unreacted  HNS.  The  principal  infra-red 
absorption  frequencies  and  the  respective  assignments  are  listed  in  Table 
1 . 

The  photographic  work  described  above  establishes  that  the  mechanical 
response  of  a  thin  layer  of  HNS  to  impact  is  characterized  by  a  large  scale 
deformation  process  prior  to  initiation  of  explosion.  As  mentioned  in  the 
Introduction,  similar  processes  are  also  observed  in  PETN,  RDX  and  HMX. 
With  PETN,  the  material  becomes  translucent  and  finally  completely 
transparent  with  rapid  flow  (300  m  s-^)  which  was  attributed  to  melting  [1], 
With  HNS  too,  high  flow  speeds  were  recorded.  A  graphical  illustration  of 
the  effect  of  mechanical  failure  of  the  sample  on  its  flow  is  provided  by 
Fig.  5  in  which  the  radius  of  the  explosive  is  plotted  as  a  function  of  time 
for  the  sequence  of  Fig.  3.  During  the  early  stages  of  compression  up  to 
100  ys  the  material  spreads  uniformly  at  10  m  s'1  and  then  the  bulk  of  it 
ceases  to  flow  till  ca.  *130  ys.  Next  there  is  a  suddent  increase  in  the 
rate  of  flow  from  zero  to  over  250  ms-1.  During  the  flow  the  layer  takes 
an  elliptical  shape  and  the  ’radius'  now  refers  to  the  major  axis  of  the 
ellipse.  Initiation  of  explosion  occurs  only  after  this  rapid  plastic 
flow. 

As  evident  from  impact  photographs  and  Fig.  14(a),  explosion  of  HNS 
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2.  Experimental 


Experiments  have  been  carried  out  using  the  drop-weight  machines 
illustrated  in  figures  1  and  2.  The  samples  used  are  5  to  7mm  in  diameter 
and  1mm  high  allowing  strain  rates  of  10^  to  1()3  s-1  to  be  obtained, 
depending  on  the  drop-height  of  the  weight.  In  both  machines  the  weight  is 
guided  by  three  cylindrical  rods.  Figure  1  illustrates  the  system  which  is 
used  to  obtain  stress-strain  curves.  Strain  gauges  attached  to  ri  allow  the 
force-time  history  to  be  monitored.  The  area  of  deforming  sample  is 
determined  by  ro  and  r^  and  the  sample  size  is  chosen  to  be,  initially, 
slightly  greater  than  this  area.  The  sample  area  under  load  is  thus 
constant  and  the  need  for  assumptions  about  volume  conservation,  which  are 
not  always  valid  for  polymers  or  porous  materials,  is  eliminated. 
Knowledge  of  the  mass  of  the  weight  combined  with  a  double  integration  of 
the  force-time  data  produces  force-displacement  information  which  is  then 
converted  to  a  stress-strain  curve.  The  system  illustrated  in  figure  2  is 
unusual  in  that  the  anvils  are  made  of  toughened  glass  and  a  light  path 
through  the  system  allows  the  sample  to  be  photographed  by  a  high  speed 
camera  at  framing  rates  up  to  2  x  1C)5  f.p.s.  during  the  course  of  the 
impact.  The  sample  strain-rates  in  both  machines  are  the  same  so  that 
direct  comparisons  between  the  photographic  records  and  the  stress-strain 
data  can  be  carried  out. 

The  third  system  used  in  this  research  was  a  miniaturised  and  modified 
Hopkinson  bar  apparatus  (as  described  by  Gorham,  1979).  This  enabled 
measurements  to  be  made,  and  stress-strain  curves  determined,  at  strain 
rates  up  to  5  x  10^  s'"1  .  The  samples  used  in  these  experiments  are  tom 
diameter  and  1mm  thick,  and  the  pressure  bar  has  a  diameter  of  3mm. 


ANVIL 


Fig.  1  Instrumented  drop-weight: 

W,  drop-weight;  H,  ri ,  r2,  r3, 
hardened  steel  rollers.  The  sample 
is  between  rollers  r2  and  r3. 


Fig.  2  Experimental  arrangement  at 
instant  of  impact:  W,  drop-weight 
(5.5  kg);  G,  glass  blocks;  M,  mirror 
P,  prism;  S,  sample. 


Finally,  the  tensile  strengths  and  rupture  strains  of  a  variety  of  PBX 
compositions  were  studied,  at  strain  rates  of  ca.  1CT11  s-1,  using  the 
’Brazilian'  test  geometry  and  laser  speckle  photography.  The  Brazilian 
test  is  illustrated  in  figure  3>  The  compression  induces  tensile  stresses 
normal  to  the  leading  axis  which  are  sensibly  constant  over  a  region  about 
the  centre  of  the  specimen.  The  tensile  strength  of  the,  material  is  then 
calculated  from 

ot  -  2P/TtDt  1 

in  which  P  is  the  failure  load,  D  and  t  the  diameter  and  thickness  of  th^ 
specincn  respectively.  The  validity  of  equation  1  is  based  on  the 
assumptions  that  failure  occurs  at  the  point  of  maximum  tensile  stress 


(i.e.,  at  the  centre),  and  that  the  compressive  stress  has  no  influence  on 
the  failure  (Fairhursl,  1 96M ) .  In  practice,  the  use  of  plane  anvils  can 
produce  very  higli  contact  stresses  at  the  loading  points,  and  lead  to  the 
collapse  of  the  contact  edge.  Awaji  and  Sata  (1979)  have  shown  that  by 
using  curved  anvils,  collapse  of  the  edge  can  be  avoided,  and  shear  stresses 
under  the  points  of  loading  may  be  substantially  reduced.  It  has  been  shown 
by  Awaji  and  Sata  that  if  the  ratio  of  the  contact  half-width,  b,  to  the 
disc  radius,  0,  is  larger  than  approximately  0.27,  the  maximum  principle 
stresses  near  the  contact  area  are  compressive.  Then  the  stresses  near  the 
centre  of  the  disc  control  the  fracture,  and  cracking  occurs  at  the  centre, 
at  a  maximum  tensile  stress  given  by 

ot*  =  (1  -  (Vr)2}  o t  •  •  •  •  2 

The  use  of  strain-gauges  to  measure  the  tensile  strain  to  failure  of  the  PBX 
compositions  have  several  disadvantages.  These  materials  are  relatively 
soft  and  the  gauge  may  reinforce  the  material  at  the  point  where  the  strain 
is  to  be  measured.  Furthermore ,  once  installed,  it  can  only  provide 
information  at  one  point.  An  alternative  approach  is  to  use  an  optical 
technique  which  measures  in-plane  surface  displacements.  Laser  speckle 
photography  is  a  simple  remote-sensing  technique,  in  which  in-plane 
displacements  can  be  determined,  with  a  high  sensitivity  Chiang  (1982). 
For  displacement  and  strain  measurements,  a  double  exposure  is 
photographically  recorded  before  and  after  deformation.  The  negative  then 
consists  of  two  identical  speckle  patterns  that  have  been  translated  with 
respect  to  eaeli  other.  Displacement  information  can  be  extracted  in  a 
point-wise  manner',  by  allowing  an  unexpanded  laser  beam  to  pass  through  the 
negative.  If  at  the  point  illuminated  by  the  beam,  the  'speckles'  have  been 
displaced  between  exposures,  the  resulting  displacment  vector  is  manifested 
by  a  Young's  fringe  pattern  observed  on  a  screen  placed  at  the  far-field 
Fraunhofer  diffraction  plane.  The  displacement  vector  is  perpendicular  to 
the  fringe  orientation  and  the  magnitude  inversl”  proportional  to  the 
fringe  spacing.  Figure  4  shows  a  speckle  pattern  being  processed . 


Fig.  3  Loading  arragement  used  in 
Brazilian  test. 

3.  Results 


laser,  L,  passes  through  the  speckle 
pattern,  P,  and  forms  a  Young's 
fringe  pattern  at  ttie  screen  S. 


3.1  pure  Explosive 

Figure  5  shows  an  impact  on  PE'IIl  taken  using  the  system  of  figure  2.  At 
firs',  the  so;.  p) e  is  compressed  anti  compacted,  and  its  radial  expansion  is  at. 
a  similar  velocity  to  that  of  impact  (i.e.  a  few  m  s”M.  This  stage  is  not 
illustrated.  Eventually  the  strength  of  the  layer  is  exceeded  and  the 
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1 1  **'•  »»>(i  i  '"cites  at.  a  number  of  Fig.  7:  (a)  Impact  on  a  disc  of  iY. 

oites.  in  fr.ai  rs  7,  ft  and  9  9  The  sample  deforms  by  bulk  plastic 

is  a  flaw  in  tin,-  glass  and  1  is  flow.  Frame  tines  ps :  a.  0;  b,  76; 

an  ignition  site.  San  pie  mass  is  c,  266.  Field  of  view  20  mm. 

1  up,,  inter-frame  time  5.5  ps.  (b)  Impact  on  a  disc  of  PS.  The 

sample  fails  catastrophically  (c). 
Frame  times  ps:  a,  0;  b,  21;  c,  28. 
Field  of  view  2.5  mm. 

material  deform;-,  plastically.  It  'is  important  to  realize  that,  in  thir 
loading  geom-fry  the  stress  for-  yield  py  is  higher  than  the  uniaxial  y  i  -  1  < 
stress  of  the  material  oy .  This  is  due  to  the  constraining  effect  of  thr 
anvils.  Application  of  the  Mises  yield  criterion  to  a  thin  layer  of 
thickness  h  and  diameter  d  between  two  rigid  cylindrical  punches  gives  the 
following  relation  (Sehroeder  and  Webster  1 9 JI9 ) 

Py  -  0 y  { 1  4  (d/3/3h)l  3 

Nrto  that.,  the  smaller  the  thickness,  or  the  greater  the  diameter,  the 
higher  the  value  for  py .  Once  plastic  flow  results,  the  radial  expansion  of 
til'-  layer  can  |x>  rapid  (more  than  l  no  ni  s’"  M .  In  explosive  samples  without 
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layer  also  goes  transparent,  but  this  is  not  due  to  fusion  but  to  sintering 
of  particles.  Evidence  for  this  is  that  a  similar  transparency  can  be 
achieved  in  AP  by  static  loading,  the  transparency  during  impact  occurs 
before  the  plastic  flow  and  is  partially  lost  during  the  plastic 
deformation  stage,  and  that  once  ignition  occurs  fractures  can  be  detected 
in  the  material  ahead  of  the  burning  front  (see  figure  4,  Heavens  and  Field 
1974).  The  observation  that  rapid  deformation  is  a  necessary  prerequisite 
for  ignition  supports  the  many  observations  of  Kholevo  (1946)  and  Afanas'ev 
and  Bobolev  (1971)  that  the  prevention  of  flow  in  an  explosive  sample  makes 
.  the  sample  less  sensitive.  Further  support  for  this  view  is  provided  by  the 
stress-time  traces  obtained  with  the  instrumented  drop-weight  machine. 
Results  for  impacts  on  a  wide  range -of  explosive  and  inert  layers  have  been 
presented  before  (Heavens  and  Field  1 97^) •  A  typical  result  is  reproduced 
in  figure  6.  The  upper  trace  represents  the  stress-time  relation  recorded 
by  the  gauges  during  an  impact  on  a  layer  of  RDX,  and  the  lower  trace 
detects  the  instant  of  explosion  by  recording  the  light  output  with  a 
photocell.  This  result  shows  clearly  the  very  close  relation  between 
mechanical  failure  of  the  layer  and  the  instant  of  ignition,  the  ignition 
■alWays  occurring  during  or  after  the  failure. 


TABLE  1 


Compound 

Melting 

point/K 

field 

stress/MPa 

H50/m 

Rapid 

flow 

Ignition 

TNT 

354 

34a 

1.48b 

Noc 

No 

Picric  acid 

395 

52 

0.73 

No 

No 

Tetryl 

404 

53 

0.37 

No 

Yesd 

Nitrocellulose 

408 

- 

• 

0.50 

No 

Yesd 

PETN 

414 

60 

0.13 

Yes 

Yes 

RDX 

474 

82 

0.28 

Yes 

Yes 

AP 

>493 

- 

- 

Yes 

Yes 

HMX 

553 

125 

0.32 

Yes 

Yes 

HNS 

588 

-140 

0.54 

Yes 

Yes 

a. '  All  the  yield  stress  values  have  been  taken  from  Afanas'ev  and  Bobolev 

(1971)  except  for  HNS  (estimated  by  Krishna  Mohan  and  Field  (1984). 

b.  Data  from  Lawrence  Livermore  Explosives  Handbook,  UCRL-52997  (1981). 

c.  The  experimental  observations  for  all  compounds  excluding  tetryl  and  HNS 
are  from  Heavens  and  Field  (1974). 

d.  Reaction  fails  to  propagate. 

Recently  we  have  shown  that  the  higher  the  yield  stress  of  a  pure  explosive 
the  higher  its  plastic  flow  velocity  in  the  drop-weight  test.  Thus  although 
with  a  strong  sample  the  onset  of  irreversible  deformation  is  delayed  it 
eventually  develops  more  rapidly  when  the  stored  energy  is  released.  This 
means  that  the  likelihood  of  ignition  during  the  conversion  of  mechanical 
energy  is  greater.  Table  1  illustrates  these  results  and  explains  why  a 
relatively  soft,  low  melting  point  explosive,  such  as  trinitrotoluene  (TNT) 
does  not  ignite  so  readily  in  the  drop-weight  test.  A  fuller  discussion  can 
be  found  in  Krishna  Mohan  and  Field  (1984). 

3.2  Explosives  with  Polymer  Additives 

Certain  polymers,  when  mixed  with  explosives,  cause  them  to  have  a  higher 
sensitiveness  to  Impact  (Swallowe  and  Field  1981,  1982).  The  evidence 
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suggests  it  is  not  due  to  chemical  effects  but  is  relatd  to  the  polymers 
thermal  properties  and  mode  of  deformation  which  both  affect  the  'hot-spot' 
formation.  Polymers  that  sensitize  are  those  that  are  relatively  tough, 
but  that  when  they  fail  do  so  catastrophically  with  the  energy  dissipation 
either  in  adiabatic  shear  bands  or  at  crack  tips.  Such  polymers  include 
polycarbonate  (PC),  polysulphone  (PS)  and  polyester  (PE).  These  materials 
all  have  low  thermal  conductivities  and  latent  heats  of  fusion  which  allow 
'hot-spot'  temperatures  to  reach  much  greater  values  than  that  of  the 
softening  point.  Conversely,  polymers  such  as  polypropylene  (PP),  which 
deform  in  a  bulk  fashion  and  which  have  relatively  high  conductivities  and 
latent  heats  have  'hot-spot'  temperatures  'quenched'  to  the  melting-point 
value.  Figure  7  illustrates  the  different  deformation  modes  of  PP  and  PC 
during  impact  in  the  drop-weight  machine.  Figure  8  illustrates  their 
stress/strain  curves,  found  using  the  apparatus  shown  in  figure  1.  The 
almost  vertical  drop  in  the  curve  for  PC  is  due  to  the  catastrophic  failure 
illustrated  in  figure  7b. 


Fig.  8  Stress-strain  curves  for  Fig.  9  Stress-strain  curves  for 

samples  of  PP  and  PC.  Catastrophic  compositions  A,  B,  C,  D,  at  a  strain 

failure  of  the  PC  sample  occurs  at  a  rate  of  2.5  x  10^  s"1. 
strain  of  ca .  1.8. 

3.3  PBX's 

Table  2  gives  information  on  the  PBX's  studied.  They  were  all  based  on  HMX 
but  with  three  different  binder  polymers.  The  amount  of  binder  and  the 
crystal  size  distributions  of  the  HMX  were  varied.  Included  in  Table  2  are 
the  results  on  the  strain  to  failure  and  the  failure  stress  obtained  using 
the  speckle  technique.  Various  general  points  can  be  made  (i)  micronized 
samples  are  stronger  (ii)  the  more  the  polymer  the  greater  the  strain  to 
failure  (iii)  PE  gives  the  most  brittle  behaviour  (lowest  strain  to 
failure)  followed  by  Viton  and  PU.  (iv)  composition  C  had  the  highest 
combination  of  strength  and  strain  to  failure.  Figure  9  and  Table  3 
summarize  results  on  the  PU  ‘compositions  from  the  high  rate  of  strain 
experiments  (c  -  2.5  x  10^  s_1).  Note  that  the  sample  with  micronized  HMX 
exhibit  the  highest  dynamic  strengths,  with  C  again  proving  the  strongest. 
Table  3  also  contains  F  of  I  and  labset  figures.  The  micronized  samples 
(A,C)  had  smaller  labset  (explosiveness)  values.  This  was  also  found  with 
the  other  binder  materials.  Figures  10  and  11  show  photographic  records  of 
Impacts  on  compositions  B  and  C.  The  micronized  PBX  (composition  C)  burns 
smoothly  and  rapidly  but  leaves  a  large  amount  of  residue.  The  sample  of 
composition  B  shows  extensive  break  up.  The  production  of  large  internal 
surfaces  and  the  rapid  penetration  of  burning  into  them  leads  to  a  more 
violent  and  complete  explosion  (less  residue  in  these  tests).  This 
probably  explains  the  higher  labset  value. 
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TABI.E  2 


Code 

Binder 

%  binder 

HMX 

%  strain 

to 

Failure  stress/ 

by  weight 

failure 

MPa. b/U  =  0.27 

A 

2 

M 

0.95 

1.38 

B 

Pol yurethanc 

2 

A/M 

0.50 

1  .22 

C 

(PU) 

5 

M 

0.98 

3.20 

D 

5 

A/M 

1 .2*1 

0.99 

E 

2 

M 

ca.  0.06 

3-11 

F 

Polyethylene 

2 

A/M 

0.09 

0.67 

G 

(PE) 

5 

M 

internal 

shear 

11 

5 

A/M 

0.18 

0.56 

"i 


J 

i 
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TABLE  3 


Code 

Dynamic  strength  (2.5  xIO^  s“1)/MPa 

F  of  Ia 

Labsetb 

A 

125 

62 

185 

B 

80 

64 

812 

C 

145 

55 

58 

D 

85 

49 

232 

a.  Figure  of  Insensitiveness,  based  on  a  drop-height  for  RDX  of  80.  The 
greater  the  figure  the  less  the  sensitiveness  (AWRE  data). 

b.  Labset  is  the  AWRE  test  for  explosiveness.  The  greater  the  figure  the 
greater  the  explosive  response  (see  Coley,  1976). 

4.  Conclusion 

The  combination  of  techniques  used  in  this  study  has  provided  valuable 
information  on  the  mechanical  and  ignition  behaviour  of  explosives.  The 

techniques  themselves  have  application  to  deformation  studies  of  materials 
in  general . 
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ABSTRACT 

This  paper  describes  a  study  of  the  decomposition  of  pentaerythritol 
tetranitrate  (PETN)  using  a  high-resolution  time-of-f light  mass 
spectrometer.  The  decomposition  was  induced  by  fracturing  single  crystals 
and  by  laser  irradiation.  In  the  fracture  experiments,  the  energy  input  was 
varied  from  the  lowest  level  necessary  to  produce  smooth  cleavage  surfaces 
to  high  energy  loading  which  produced  rough  concholdal  fracture  surfaces. 
In  the  laser  experiments,  a  ruby  laser  was  used  in  both  normal  and 
Q-switched  modes,  and  again  the  energy  input  wa3  varied.  For  all  the 
various  experiments,  the  reaction  products  were  analysed  and  reaction 
schemes  are  proposed.  It  is  shown  that  low  energy  fracture  causes 
decomposition  which  follows  the  same  reaction  pathway  as  that  induced 
thermally,  with  initial  failure  at  the  RO-NO2  bond.  However,  high  energy 
fracture  results  in  the  breaking  of  the  C-C  bonds.  Two  reaction  pathways 
were  observed  with  the  laser  irradiation.  The  first  is  the  normal  thermal 
process,  but  evidence  was  also  found  for  failure  at  the  R-ONO2  bond.  The 
reaction  continued  for  several  milliseconds  after  the  end  of  the  laser 
pulse  suggesting  a  'partial'  ignition  of  the  explosive.  In  other 
experiments,  the  conditions  for  laser  initiation  of  PETN  in  vacuum  were 
investigated.  Explosion  occurred  when  Q-switched  pulses  of  1  J  energy  were 
applied  to  a  molten  layer  of  PETN. 
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1  INTRODUCTION 

Pentaerythritol  tetranitrate  (PETN)  is  a  widely  used  secondary 
explosive.  Its  relatively  high  chemical  stability  allows  it  to  oe 
manufactured,  handled  and  stored  in  comparative  safety.  Its  thermal 
decomposition  has  been  the  subject  of  intensive  study,  and  results  on  the 
energies  involved  and  the  reaction  pathway  are  summarised  below  in 
section  2.  Such  decomposition  has  relevance  to  the  long-time  storage  of 
PETN  and  to  situations  involving  relatively  slow  thermal  inputs.  However, 
there  is  no  certainty  that  the  same  reaction  pathways  will  operate  when  fast 
decomposition  at  transient  'hot  spots*  takes  place. 

It  is  now  generally  accepted  that  explosives  almost  Invariably  ignite 
by  thermal  processes  (Bowden  and  Yoffe  1952,  1958).  The  energy  input  may  be 
directly  thermal  or  may  involve  the  degradation  of  mechanical  energy  into 
local  'hot  spots'  (see  Field,  Swallowe  and  Heavens  1983,  for  a  recent 
review).  Whether  or  not  fast  reaction  develops  from  the  'hot  spot'  then 
depends  on  the  balance  of  heat  produced  by  chemical  reaction  and  that  lost 
by  dissipative  processes 

The  object  of  the  present  research  was  to  investigate  the  decomposition 
pathways  fo;  PETN  following  the  production  of  transient  'hot  spots' ,  and  to 
compare  the  results  with  those  obtained  from  slow  thermal  decomposition 
studies.  In  one  set  of  experiments,  crystals  were  fractured  in  a  variety  of 
ways;  in  this  situation  the  'hot  spot’  forms  at  the  tip  of  the  running  crack 
(see  Field  et  al  (1982)  for  a  discussion).  Direct  measurements  of  crack  tip 
temperatures ,  which  may  reach  several  hundred  K  have  been  made  by  Doll 
(1973),  Fuller,  Fox  and  Field  (1975)  and  Weichert  and  Schonert  ( 1 97 if) .  In 
other  experiments,  a  ruby  laser  was  used  to  produce  a  transient  'hot  spot' 
in  the  material.  The  experiments  were  performed  in  an  ultra-high  vacuum 
system  with  a  time-of-f light  (TOF)  mass  spectrometer  used  to  detect  and 
analyse  the  transient  species  and  the  decomposition  products.  A 


preliminary  account  of  the  fracture  experiments  was  given  in  Hauser,  Field 
and  Krishna  Mohan  (1983). 

Finally,  the  conditions  under  which  PETN  initiates  in  vacuum  when 
irradiated  by  a  ruby  laser  pulse  were  investigated. 

2  THERMAL  DECOMPOSITION 

Andrews  and  Kaidymov  (I960  found  an  activation  energy  E  of  215  to 
263  kJ  mol-1  for  the  decomposition  of  solid  PETN  but  their  result  was  based 
on  a  limited  amount  of  data.  Values  from  other  workers  using  a  variety  of 
methods  fall  reasonably  close  together.  Rogers  and  Morris  (1966)  using 
differential  scanning  calorimetry  (DSC)  found  a  value  of  198  kJ  mol-1. 
Soria-Ruiz  (1969)  obtained  176  kJ  mol-1  by  monitoring  the  peak  m/e  28  using 
a  quadrapole  mass  spectrometer.  Ng,  Field  and  Hauser  (1976),  measuring  the 
total  gas  evolution,  found  a  value  of  192  kJ  mol-1.  Vapor  phase 
decomposition  was  studied  by  Robertson  (19*18)  who  found  an  activation 
energy  of  197  kJ  mol-1  while  Maycok  and  Verneker  (1970)  using  a 
thermobarogravimetri c  technique  obtained  1 88  kJ  mol-1.  Decomposition  is 
complex  and  takes  place  concurrently  with  sublimation.  Decomposition 
products  at  *183  K  were  analysed  chemically  by  Rideal  and  Robertson  (19*18) 
who  found 

CjHgNuO^  -  0.52N02  ♦  2.  UNO  ♦  0.»J2N20  +  0.28C02  +  0.93C0  +  0.09H2  +  water  + 
formaldehyde  ♦  residue  (1) 

This  was  supported  by  the  time-of-f light  mass  spectrometry  study  of  Ng  et  al 
(1976)  who  found  that  the  product  distribution  is  temperature  sensitive. 
The  activation  energy  for  decomposition  is  close  to  the  bond  dissociation 
energy  of  ca.  173  kJ  mol-1  for  R0-N02  in  mononitrate  and  nitrite  esters 
estimated  by  Levy  (195*1).  A  mechanism  was  proposed  by  Ng  et  al  (1976)  for 
the  thermal  decomposition  of  PETN  via  the  rupture  of  the  R0-N02  bond 
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followed  by  the  elimination  of  a  neutral  formaldehyde  molecule  as  given 
below: 

Normal  Thermal  Decomposition  Pathway 


C(CH20N02)4 

-► 

C(CH^m02)^CH^)  +  no2 

2 

C  (CH2ONO  2)3.011 2*3 

.C(CH20N02)3  +  HCH0 

3 

.C(CH20N02)3 

+ 

2CH50N02  +  2C0  +  NO 

1» 

CH3ONO2 

ch3o  +  N02 

5 

CH3O  +  no2 

-*• 

HCH0  +  H0N0 

6 

2H0N0 

-► 

♦  NO  +  N02 

7 

CH3O  +  HCH0 

CH^H  ♦  HCO 

8 

HCO  +  NO 

CO  +  HN0 

9 

2HN0 

♦  N^ 

10 

3  EXPERIMENTAL 

PETN  in  powder  form  was  obtained  from  PERME,  Waltham  Abbey.  Large 
single  crystals  (typically  20  x  5  x  3  mm3)  were  grown  from  acetone  solution. 
The  crystals  were  prismatic  with  elongated  hexagonal  (100)  planes.  They 
were  mounted  on  a  movable  platform  inside  a  six-way  reaction  chamber  of  an 
ultra-high  vacuum  (UHV)  system  connected  to  a  Bendix  RGA-1A  time-of-flight 
mass  spectrometer.  A  view  of  the  experimental  arrangement  is  given  in 
figure  1.  The  sample  was  within  40  mm  of  the  ion  source  which  was  operated 
at  70  V. 

A  trace  offset  controller  was  used  to  present  a  step-wise  display  of  the 
mass  spectra  on  a  Tetronlx  oscilloscope  (model  A7603).  It  is  capable  of 
displaying  a  maximum  number  of  64  spectra  in  8  steps  on  a  single 
oscillogram.  The  sweep  time  for  each  spectrum  is  33  ps.  It  also  has  a 
horizontal  offset  component  to  displace  the  steps  sidewise  relative  to  each 
other  so  as  to  avoid  overlapping  of  high  rising  peaks  between  adjacent 
steps.  A  built-in  time  delay  unit  on  the  oscilloscope  trace  offset 


controller  Is  adjustable  between  HpS  and  1  ms.  An  additional  time  delay 
circuit  was  designed  (Hauser  1977)  to  extend  and  vary  the  delay  time  between 
HO  ms  and  3  s.  A  high-writing  speed  camera,  C27,  using  polaroid  film  was 
used  to  record  the  oscillograms.  An  oscilloscope  bright-up  unit  was  also 
used  to  facilitate  easy  data  acquisition, 
a)  Fracture  Induced  Decomposition 

A  chisel  was  mounted  on  a  metal  bellows  assembly  situated  on  the  top 
compartment  of  a  six-way  reaction  vessel.  Its  shaft  had  a  hollow  section  to 
accommodate  a  small  explosive  fuse.  The  platform  on  which  the  sample  was 
placed  could  be  raised  towards  the  stationary  chisel  edge  by  means  of  a 
linear  motion  device  (LMD)  mounted  in  the  compartment  opposite  to  the 
chisel.  When  the  crystal  surface  was  just  touching  the  chisel  edge,  two 
modes  of  initiating  a  fracture  were  then  carried  out. 

i )  Less  energetic  fracture 

The  chisel  edge  was  aligned  parallel  to  a  (110)  cleavage  plane  of  the 
crystal  via  a  rotating  flange  to  which  the  LMD  was  attached  (figure  2a). 
This  could  be  achieved  by  rotating  the  entire  LMD  assembly  until  the  edge  of 
the  chisel  wa3  parallel  to  the  long  edge  of  the  cyrstal.  The  crystal  was 
loaded  against  the  chisel  by  raising  the  platform  until  the  crystal 
cleaved.  This  process  was  carried  out  slowly  so  that  the  minimum  cleavage 
force  was  used.  A  spectrum  was  taken  as  soon  after  cleavage  as  possible. 
The  estimated  time  delay  between  the  event  and  the  recording  was  ca.  1  s. 

i i )  Energetic  fracture 

In  this  situation,  the  crystal  was  not  orientated  preferentially  and 
excess  energy  was  rapidly  fed  into  the  chisel  by  firing  an  explosive  fuse 
placed  in  its  shaft  (see  figure  2b).  In  these  experiments,  synchronization 
of  the  firing  of  the  detonator  with  the  mass  spectrometer  recording  was 
achieved  and  chosen  delay  times  were  possible. 


b)  Laser-induced  Decomposition 

Laser  irradiation  of  a  sample  in  the  chamber  was  achieved  by  using  the 
transparent  windows.  The  Bradley  ruby  laser  (A  =  69^.3  nm)  was  first 
aligned  with  the  help  of  a  small  He-Ne  continuous  laser  which  was  mounted  on 
the  same  optical  bench.  A  lens  was  used  to  focus  the  ruby  beam  onto  the 
crystal  surface  to  an  area  of  ca.  2.0  ±  0.5  mm  in  diameter.  An  ion-repeller 
grid  operating  at  positive  voltage  was  situated  between  the  sample  and  the 
ion  source.  When  operated  in  Q-switched  mode,  the  ruby  laser  has  a  peak 
power  of  ca.  200  MW  and  a  pulse  duration  of  ca.  20  ns.  In  the  decomposition 
experiments,  energy  inputs  into  the  crystals  were  in  the  range  of  0.1  to 
1  J.  The  valve  leading  to  a  Ferranti  ion  pump  was  adjusted  such  that  the 
oressure  inside  the  system  did  not  exceed  10~5  torr  at  maximum  gas 
evolution . 

c )  Laser  Initiation 

In  a  separate  series  of  experiments,  the  conditions  for  the  laser 
initiation  of  PETN  in  the  vacuum  system  were  investigated.  Various 
experimental  arrangements  were  used.  In  experiments  with  solid  PETN, 
single  crystals  with  their  back  surface  coated  with  a  thin  film  (ca.  100  nm) 
of  aluminium  were  used.  The  crystals,  of  millimeter  dimensions,  were 
either  mounted  freely  (unconfined)  on  a  support,  or  embedded  in  epoxy  resin 
to  study  the  effect  of  confinement.  Molten  PETN  was  obtained  using  the 
arrangement  illustrated  in  figure  3.  The  aluminium  block  was  in  contact 
with  -i  hot  plate.  When  the  temperature  was  raised  above  the  melting  point 
(  ■-I'.U  Y  )  ,  the  melt  was  trapped  between  the  microscope  slide  and  the  aluminium 
block.  In  an  alternative  apparatus  a  molten  layer  was  produced  without  the 
confining  glass  slide. 

;i  RESULTS 

a  i  Fracture-induced  Decomposition 

The  fracture  of  a  PETN  crystal  by  an  explosively  driven  chisel  created 
new  surfaces  together  with  the  rapid  evolution  of  gases.  The  fracture 
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surfaces  for  this  energetic  mode  of  fracture  were  invariably  rough 
(figure  8b).  Their  mass  spectra  showed  distinct  features  different  from 
that  of  thermal  decomposition.  In  the  low  mass  region,  peaks  at  m/e  16(0), 
28(CO),  30(N0),  JJiKCOp.NpO)  were  produced.  However,  no  peak  at  m/e  46(N02) 
was  observed.  In  the  high  mass  regions,  the  only  prominent  peak  was  at 
m/e  60.  Figure  4  is  a  step-wise  presentation  of  the  mass  spectra  taken 
700ms  after  the  firing  of  the  detonator.  The  entire  mass  spectra  of  the 
reaction  products  are  given  in  figure  5  which  was  taken  with  a  time  delay  of 
350  us. 

In  contrast,  the  low  energy  fracture  mode,  achieved  by  the  method  of 
figure  2a,  produced  a  number  of  different  spectra  which  depended  on  the 
amount  of  energy  transmitted  to  the  crack.  Three  types  of  spectra  were 
observed . 

i)  a  large  intense  peak  at  m/e  76  plus  other  weak  ones.  An  example  is 
shown  in  figure  6. 

ii)  a  large  intense  peak  at  m/e  60  plus  others.  This  type  is  similar 
to  the  example  given  in  figure  4. 

iii)  the  simultaneous  presence  of  both  peaks  at  m/e  76  and  60  in 
addition  to  others.  An  example  is  illustrated  in  figure  7. 

The  presence  of  a  peak  at  m/e  76  is  also  accompanied  by  a  peak  at  m/e  46  in 
the  low  mass  region.  The  observed  fracture  surfaces  were  smooth  (figure  8a 
gives  an  example)  for  spectra  (i),  rough  for  spectra  (ii)  and  partly  smooth 
and  partly  rough  for  spectra  (iii).  The  smooth  surfaces  were  (110)  cleavage 
faces  and  were  obtained  when  the  chisel  was  accurately  aligned  and  the 
minimum  of  energy  used, 
b )  Laser- i nduced  Decomposition 

The  interaction  of  a  Q-switched  ruby  laser  pulse  with  a  PETN  crystal 
inside  the  IJHV  system  produced  positively  charged  species  which  were 
detected  for  a  duration  of  ca.  2  ms  after  the  firing  of  the  laser  pulse  even 
without  the  ion-3ource  of  the  spectrometer  in  operation.  A  repeller 
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-  -  x,  y  *  -  z;  x  *  -,  -  -  y,  -  -  z;  y  +  -,  x  +  ^  z  + 
l_y  1-X  z  +  1 

The  ^etrago^al  unit  cefl  which  contains  two  molecules  of  PETN  is  depicted  in 
figure  12  which  shows  projection  onto  the  (001)  plane.  The  cleavage  planes 
are  (110)  and  (110)  since  only  the  weak  van  der  Waals  forces  need  to  be 
overcome  along  these  planes.  Figure  13  shows  the  position  of  the  cleavage 
plane  (110)  relative  to  another  cleavage  plane  (110)  which  is  also  one  of 
the  surfaces  of  the  as-grown  crystal.  However,  the  rapid  fracture, 
initiated  by  the  explosively-driven  chisel,  does  not  remain  on  the  planes 
but  penetrates  and  bifurcates  into  many  non-cleavage  planes  thus  producing 
a  rough  conchoidal  surface  (figure  8b).  As  the  crack  propagates  along  the 
terrain  of  this  rough  surface,  it  has  to  overcome  and  break  off  all  chemical 
bonds  that  are  transverse  to  its  path.  However,  if  fracture  is  confined 
along  the  cleavage  plane,  then  no  such  rupture  of  chemical  bonds  takes 
place.  The  fractured  surfaces  are  then  smooth  as  illustrated  in  figure  8a. 
In  this  latter  event  the  energy  at  the  crack  tip  dissipates  in  the  form  of 
heat  causing  those  molecules  at  the  tip  to  sublime  and  subsequently 
decompose  in  the  vapour  phase.  The  mass  spectrum  is  similar  to  the  one 
obtained  from  thermal  decomposition  in  which  a  PETN  molecule  decomposes  via 
a  molecular  elimination  mechanism  (Ng  et  al  1976).  The  prominent  peak  in 
the  high  mass  region  at  m/e  76  is  attributed  to  the  species  CH20NC>2+  derived 
from  the  parent  moleule  CH3ONO2.  The  distinct  difference  in  the  mass 
spectrum  of  the  decomposition  products  and  its  correlation  with  the  mode  of 
fracturing  strongly  suggest  that  there  are  two  different  decomposition 
pathways  for  the  fracture-induced  decomposition  of  PETN.  Indeed  both  the 
reaction  pathways  can  occur  simultaneously  if  fracture  occurs  partly  along 
the  cleavage  plane  and  partly  by  bifurcating  into  many  non-cleavage  planes. 
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In  such  cases  the  mass  spectrum  shows  the  presence  of  both  the  peaks  at  m/e 
60  and  m/e  76  (see  figure  7)  in  the  high  mass  region. 


cleavage  at  the  four  central  C-C  bonds  and  that  the  fragments  .CH2ONO2  and 
.C(CH20N02)2  rapidly  degrade  to  CO  and  NO. 

In  the  mass  spectrometer ,  the  species  CH20N02+*  if  present,  should 
appear  at  m/e  76.  The  fact  that  it  is  absent  in  our  TOF  mass  spectra  for 
enegetic  fracture  indicates  that  this  radical  is  a  short-lived  species  with 
a  life-time  less  than  the  time  resolution,  33  ps,  of  the  spectrometer. 
There  are  two  degradation  pathways  for  this  radical 

i)  .CH2ONO2  -*  .CH2O  +  N02 

ii)  .CH2ONO2  +  .CH2ONO  +  0 

The  absence  of  a  peak  at  m/e  H6  (N02+)  rules  out  pathway  (i).  The  presence 
of  a  peak  at  m/e  60  (.CH20N0+)  and  m/e  30  (N0+)  (which  reflects  easy 
cleavage  of  the  alkyl  nitrite  at  R0-N0)  supports  pathway  (ii). 

The  bond  dissociation  energy  for  a  typical  (CH 3 ) 3C-CH3  molecule  is 
290  kJ  mol-1  (Benson  and  O'Neal  1970)  and  for  RO-NO2  is  173  kJ  mol-1  (Levy 
195*1,  1956).  When  some  of  the  elastic  strain  energy  dissipated  at  the  crack 
tip  goes  into  promoting  chemical  decomposition,  the  molecule  should 
decompose  via  the  rupture  of  the  weaker  bond.  Our  result  showed  that  it  is 
the  stronger  C-C  bond  that  is  broken  whereas  the  weaker  RO-NO2  bond  remains 
intact.  This  implies  that  the  energy  associated  with  the  crack  tip  is 
localised  specifically  on  the  C-C  bond  and  that  this  bond  is  traversed  by 
the  crack  in  a  time  comparable  to  that  required  for  one  vibration.  For  this 
reason  there  is  no  possibility  for  the  molecu"  to  attain  thermal 
equilibrium. 

PETN  crystallises  from  acetone  in  a  tetragonal  crystal  system  with 
space  group  P42^c  (Wyckoff  I960).  The  molecule  which  has  Si|  symmetry  has 
four  identical  CH2ONO2  groups  attached  to  the  central  carbon  which  occupies 
special  positions,  0,  0,  0  and  1/2,  1/2,  1/2.  The  other  carbon  atoms  are 
situated  at  the  general  positions  x,  y,  z;  x,  y,  z;  y,  x,  z;  y,  x,  z; 
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Unconfined  crystals  simply  lost  the  aluminium  film  during  laser 
irradiation  and  showed  only  a  small  crate’  formed  at  the  focal  point. 
Crystals  confined  in  epoxy  resin,  were  shattered  into  small  pieces  and  it  is 
probable  that  initiation  could  have  been  achieved  with  a  slightly  higher 
laser  energy. 

5  DISCUSSION 

a)  Fracture-induced  Decomposition 

When  a  fast  moving  crack  propagates  through  a  solid,  the  condition  for 
crack  growth  in  terms  of  a  Griffith-type  energy  balance  is  given  by  (Bowden, 
Fox  and  Soria-Ruiz  1968,  Fox  and  Soria-Ruiz  1970): 


d{Ee  -  [Es  +  Ek  +  Ep  +  Ec] }/dc  >  0 


where  c  is  the  crack  length  and  the  subscr  pts  for  the  energy  (E)  terms 
represent  the  contribution  of  elastic,  surface,  kinetic,  plastic  and 
chemical  energies  respectively.  Fracture  studies  on  PETN  crystal  (Kagan 
and  Chaudhri  1977)  have  shown  that  it  is  a  weak  solid  which  undergoes 
brittle  fracture  with  limited  plastic  flow.  Its  fracture  surface  energy  of 
0.11  J  m~2  is  typical  of  that  for  a  van  der  Waals  solid.  A  characteristic 
of  brittle  fracture  is  that  cracks  can  reach  very  high  velocities  of  the 
order  of  a  thousand  metres  per  second  depending  on  the  solid  (see,  for 
example,  Field  1970  for  a  review).  The  crack  tip  radius  is  of  the  order  of 
atomic  dimensions,  and  when  it  develops  it  does  so  by  a  process  of  breaking 
atomic  bonds.  In  their  work  on  fracture-induced  decomposition  in  brittle 
crystalline  solids  Bowden,  Fox  and  Soria-Ruiz  reported  the  evolution  of 
single  gaseous  species  from  metal  carbonates  and  azides  as  a  result  of 
fracturing  these  solids.  However  in  the  more  complex  PETN  crystal  there  is 
more  than  one  type  of  chemical  bonding  present  in  the  molecule.  Their 
preliminary  results  suggested  that  decomposition  takes  place  via  the 
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voltage  of  3  V  was  needed  to  stop  these  charged  species  from  entering  into 
the  ion  source  to  observe  the  mass  spectra  of  the  decomposition  products. 

With  the  ion  repeller  set  at  3  V  DC ,  the  ionisation  potential  at  70  eV 
and  the  time  delay  at  1.25  ms  the  mass  spectrum  (shown  in  figure  9)  revealed 
the  presence  of  decomposition  products  similar  to  that  obtained  from  our 
thermal  decomposition  work.  With  a  longer  time  delay  of  3.3  ms,  more  peaks 
(m/e,  60,  62,  65,  71,  7*1,  76,  9*4,  114,  116,  128,  190,  etc)  were  observed  in 
the  high  mass  region  (figure  10).  All  of  these  were  transient  with  the 
exception  of  m/e  76  which  was  particularly  persistent.  At  a  time  delay  of 
5.6  ms  many  of  those  peaks  were  still  present  but  with  diminishing  intensity 
whereas  the  peak  at  m/e  76  appears  in  every  step  (figure  11). 
c)  Laser  Initiation 

The  results  of  the  laser  initiation  experiments  are  summarised  in  table 
1.  Only  under  one  set  of  conditions  (PETN  in  liquid  phase,  confined  and 
irradiation  by  a  Q-switched  laser  pulse)  could  initiation  of  PETN  by  ruby 
laser  irradiation  in  vacuum  be  effected  for  the  energies  used.  Solid  PETN 
could  not  be  initiated  by  laser  irradiation  under  vacuum  conditions  for  the 
range  of  energies  available. 


TABLE  1 

Laser  initiation  of  PETN  in  vacuum 


PETN  phase 

Confinement 

Laser  pulse 

Explosion 

Solid 

Yes 

Non  Q-switched 

(20  J) 

No 

Solid 

Yes 

Q-switched 

(1  J) 

No 

Solid 

No 

Non  Q-switched 

(20  J) 

No 

Solid 

No 

Q-switched 

(1  J) 

No 

Liquid 

No 

Non  Q-switched 

(20  J) 

No 

Liquid 

No 

Q-switched 

(1  J) 

No 

Liquid 

Yes 

Non  Q-switched 

(20  J) 

No 

Liquid 

Yes 

Q-switched 

(1  J) 

Yes 

i 


Jkii 


Slow  thermal;  low  energy  fracture 
low  energy  laser. 

High  energy  laser. 


46 


5.  The  deformation  stress  of  highly  brittle 
explosive  crystals  from  real  contact  area 
measurements 

M.  MUNAWAR  CHAUDHRI 

Cavendish  Laboratory,  Madingley  Road,  Cambridge,  Uh' 


An  experimental  study  is  described  in  which  the  deformation  stress  of  millimetre-size 
highly  brittle  as-grown  single  crystals  of  explosive  materials  has  been  determined  using  a 
novel  experimental  technique.  The  technique  involves  the  loading  of  the  test  crystal 
against  a  transparent  glass  plate  with  in  situ  measurements  of  the  real  area  of  contact 
under  load  using  an  optical  microscope.  For  all  explosive  crystals  the  deformation  stress 
was  found  to  increase  with  load  reaching  a  plateau  at  higher  loads.  The  maximum  defor¬ 
mation  stress  of  a  crystal  was  found  to  be  considerably  lower  than  its  Vickers  diamond 
hardness  value.  A  discussion  on  the  determination  of  the  uniaxial  yield  stress  of  small 
particles  is  also  given. 


1.  Introduction 

Recently,  Chaudhri  [1]  published  a  model  for  the 
initiation  of  explosions  in  compacts  of  sensitive 
solid  explosives  when  they  are  impacted  by  hard 
needle-shaped  strikers.  According  to  this  model, 
an  explosion  results  due  to  the  high  temperatures 
generated  by  friction  at  the  contacting  surfaces  of 
the  explosive  particles  which  are  picked  up  by  the 
striker’s  tip  and  carried  into  the  explosive  compact. 
To  make  reliable  quantitative  estimates  of  the  tem¬ 
peratures  generated  at  the  rubbing  surfaces,  a 
knowledge  of  several  physical  properties  of  the 
particles  is  needed  [1].  One  such  property  of 
significant  importance  is  the  deformation  stress 
of  a  particle  pressed  against  another  with  a  given 
normal  load.  To  determine  this  parameter  direct 
measurements  of  the  real  area  of  contact  are 
required,  which  are  difficult  to  make  even  for 
ductile  metals  as  has  recently  been  pointed  out 
by  Tabor  [2j.  For  brittle  materials  the  situation 
is  further  complicated  by  cracking  which  can 
occur  at  quite  small  loads.  Here  the  results  of  a 
novel  experimental  investigation  are  reported  in 
which  the  deformation  stress  of  a  number  of 
highly  brittle  [3,  4]  as-grown  explosive  crystals 
has  been  determined  from  measurements  of  the 
real  area  of  contact  when  they  are  pressed  against 


a  hard  transparent  surface  with  a  known  normal 
load.  In  the  literature,  such  a  method  of  deter¬ 
mining  the  deformation  stress  of  brittle  crystals 
of  shapes  other  than  the  rectangular  parallelopiped 
does  not  appear  to  have  been  reported  before. 

2.  Experimental  details 
2.1.  Materials 

Single  crystals  of  lead  azotetrazole  (LAT)  of 
types  A  and  C  were  obtained  from  Dr  P.  Collins 
of  PERME  (UK);  type  A  were  of  pyramidal  shape 
about  2  mm  across,  whereas  type  C  were  needle-like 
1  to  2  mm  in  length  and  about  O.S  mm  in  diameter. 
Single  crystals  of  a-lead  azide  were  grown  in  the 
laboratory  by  a  diffusion  method  [5];  these  were 
also  needle  shaped,  0.3  to  0.5  mm  in  diameter, 
and  1  to  2  mm  in  length.  Crystals  of  0-cyclotetra- 
methyl-tetranitramine  (/3-HMX)  and  cyclotri- 
methylenetrinitramine  (RDX)  were  also  grown  in 
the  laboratory  by  slow  evaporation  of  a  saturated 
solution  in  acetone  at  70°  C.  0-HMX  single  crystals 
were  thick  platelets  with  relatively  good  edges, 
whereas  crystals  of  RDX  did  not  show  sharp 
edges  and  generally  had  rounded  corners.  Potassium 
chloride  crystals  about  1  cm  x  1  cm  x  1  cm  in  size 
were  obtained  from  BDH  (UK),  from  which  pieces 
about  1  to  2  mm  in  size  were  cleaved  for  these 
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TABLE  1  Chemical  formuli  and  some  crystallographic  and  mechanical  properties  of  the  crystals  used  in  this  work 


Material 

Chemical  formula 

Crystal  class 
and  space  group 

Primary  glide 
system 

o-lead  azide 
(o-PbN.) 

PbN, 

Orthorhombic 

Pcmn  (6) 

not  known 

Lead  azote tr azole  (LAT) 

type  A 

(C,N„,Pb]x  •  H,0 
(x  between  0  and  2) 

not  known 

not  known 

type  C 

(C,N,0Pb]  ■  Pb(OH), 

not  known* 

not  known 

3-cyclotetramethylene- 
tetranitramine  (0-HMX) 

<CH,NNO,)« 

monoclinic 

/*2,/e|8) 

not  known 

Cyclotrimethylene- 

C)H4N,(NO,), 

orthorhombic 

(0  1 0)(00  1 1 

trinitramine  (RDX) 

Pbca  (9) 

[10] 

Potassium  chloride 

KC1 

cubic 

{1  10}ul0> 

(KC1) 

Fm 3F  J 1 1 ) 

1121 

‘Melanie  A.  Pierce-Butler  |7]  has  determined  the  structure  of  LAT  type  D  having  the  same  chemical  formula  as  that  of 
the  type  C;  the  crystal  class  is  triclinic  and  space  group  PI. 


experiments.  Table  I  lists  the  chemical  formulae  of 
the  various  explosives  along  with  their  crystal 
structures  and  slip  systems  (if  known). 

2.2.  The  principle  of  the  experimental 
technique 

The  principle  of  the  technique  can  be  illustrated 
with  the  help  of  Fig.  1.  In  Fig.  la  a  light  beam  is 
reflected  from  a  partially  silvered  surface  S  and 
hits  a  polished  glass  plate  normally.  If  the  refrac¬ 
tive  indices  of  the  glass  and  air  are  nt  and 
respectively,  then  the  intensity  /,  of  the  light 
reflected  from  surface  A,  for  unit  incident  inten¬ 
sity,  is: 


Substituting  n  t  =  1.52  and  rt3  -  l,  we  have  l\  — 
0.04. 

A  crystal  of  the  test  material  is  then  brought 


into  perfect  contact  at  surface  A  as  shown  in 
Fig.  lb.  The  intensity  /2  of  the  light  reflected  from 
the  glass/crystal  interface  is  then  given  by: 


where  n3  is  the  refractive  index  of  the  crystal. 
For  the  crystals  used  in  this  work  n3  is  in  the 
range  1 .5  to  2.6.  Taking  a  value  of  n3  of  1.6  (say), 
/j  =  0.0066,  which  is  considerably  smaller  than 
/ 1.  The  same  is  true  for  other  values  of  n3  in  the 
given  range.  Thus,  the  reduction  in  the  intensity 
of  the  reflected  light  gives  is  a  convenient  method 
of  measuring  the  area  of  eal  contact  between  a 
non-metallic  crystal  and  a  tiansparent  glass  plate. 

2.3.  The  compression  apparatus 
A  schematic  diagram  of  the  arrangement  used  for 
compressing  crystals  against  an  optically  flat  glass 
surface  and  simultaneously  measuring  the  area  of 


Figure  I  Diagram  showing  the  principle 
of  the  experimental  technique  used  for 
measuring  the  real  area  of  contact. 
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Figure  2  Schematic  diagram  of  the  arrangement  used  for 
measuring  the  real  area  of  contact  of  a  crystal  loaded 
against  a  smooth  glass  plate. 

real  contact  between  the  crystal  and  the  glass 
surface  is  shown  in  Fig.  2.  The  crystal  to  be 
examined  is  mounted  on  a  steel  holder  with  an 
epoxy  “araldite”.  Before  the  epoxy  hardens,  the 
crystal  is  carefully  positioned  so  that  as  far  as 
possible  a  sharp  corner  of  the  crystal  is  along  the 
axis  of  the  holder,  or  a  sharp  edge  is  normal  to  it. 

The  steel  holder  itself  is  screwed  into  a  sensitive 
load  cell  capable  of  measuring  loads  from  less  than 
a  gramme  to  several  kilogrammes.  The  load  cell  is 
rigidly  mounted  on  the  crosshead  of  a  mechanical 
testing  machine  (Instron  No.  1122);  in  all  exper¬ 
iments  the  crosshead  speed  was  0.05  mm  min'1. 
The  glass  plate  is  placed  on  a  thick  steel  support 
plate  and  the  contact  area  between  the  crystal  and 
the  glass  plate  is  viewed  in  reflection  using  an 
optical  microscope  focused  on  the  glass  surface, 
as  shown  in  the  figure.  When  desired,  photographs 
of  the  contact  area  are  taken,  the  magnification 
on  the  camera  film  being  x  60. 

The  experimental  procedure  was  as  follows. 
First  the  crystal  was  brought  to  within  5  to  10/rm 
of  the  glass  surface  and  a  photograph  taken.  At 
this  stage  the  field  of  view  was  uniformly  illumi¬ 
nated.  Then  the  crystal  was  brought  into  contact 
with  the  glass  surface,  with  the  normal  load  being 
a  gramme  or  less.  It  was  observed  that  wherever 
the  crystal  made  real  contact  with  the  glass  sur¬ 
face,  the  intensity  of  illumination  decreased 
markedly.  The  load  was  then  further  increased 
by  a  small  amount  and  another  photograph  taken. 
This  procedure  was  continued  for  several  loads 
and  a  corresponding  sequence  of  photographs  was 
made. 

The  maximum  load  applied  on  a  crystal 
depended  upon  its  size  and  the  type  so  that  only 


a  minimum  amount  of  cracking  occurred.  Photo¬ 
graphs  were  also  taken  during  the  unloading  of  the 
crystal  and  it  was  found  that  in  all  cases  plastic 
deformation  of  the  crystal  had  occurred  over  the 
entire  region  in  real  contact  with  the  glass  plate. 
The  deformed  region  was  flat  and  optically  smooth 
and  when  the  crystal  was  only  just  removed  from 
the  glass  plate  the  flat  gave  a  mirror-like  reflection. 

Silhouette  optical  micrographs  of  crystals 
before  and  after  deformation  were  also  taken. 
Although  the  deformed  regions  were  relatively 
small,  in  some  cases  the  deformation  could  be 
clearly  seen  in  the  migrographs. 

The  measurements  of  the  real  area  of  contact 
were  made  from  enlarged  prints  using  a  transparent 
millimetre  grid.  In  this  manner,  a  real  contact  area 
as  small  as  1  /urn  x  1  #nn  could  be  measured. 

3.  Results 

The  results  from  different  crystals  are  given  below 
separately. 

3.1.a-PbN6 

Two  single  crystals  of  this  material  were  investi¬ 
gated.  Fig.  3a  shows  the  profile  of  a  crystal  before 
it  was  deformed.  It  was  loaded  on  an  edge  (marked 
with  arrow)  along  the  c-axis,  and  the  photographs 
of  real  contact  are  shown  in  Fig.  3b.  The  load  on 
the  crystal  in  grammes  is  shown  alongside  each 
frame.  In  the  uppermost  frame,  the  crystal  has  not 
touched  the  glass  plate  and  is  only  a  few  /an  from 
the  surface.  Note  that  the  entire  field  is  uniformly 
illuminated  (the  black  circular  spot  in  the  centre 
of  frames  1  and  2  is  due  to  a  dust  particle).  In 
frame  2,  the  load  has  been  increased  to  10.5  g 
and  a  dark  band  has  appeared  in  the  centre.  The 
crystal  is  in  real  contact  with  the  glass  surface  at 
all  points  in  the  band.  In  the  next  frame,  the 
load  has  been  further  increased  to  20  g,  and  note 
that  the  area  of  the  dark  band  has  also  increased. 
The  crystal  was  then  unloaded,  first  to  4g  and 
then  to  zero.  At  4  g  the  area  of  the  dark  region  in 
the  band  has  decreased  as  some  of  it  has  become 
brighter  than  the  background.  At  zero  load,  the 
entire  band  has  changed  into  a  bright  one.  These 
observations  strongly  indicate  that  the  crystal  is 
no  longer  in  contact  with  the  glass  plate  and  that 
the  edge  of  the  crystal  had  plastically  deformed  to 
a  smooth  flat  surface. 

The  crystal  was  then  again  loaded  and  the  real 
areas  of  contact  corresponding  to  loads  of  24,  32 
and  41  g  are  also  shown  in  the  figure.  As  the  load 
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Figure  3  (a)  Profile  of  an  a-lead  azide  single  crystal  which 
was  loaded  along  the  c-axis;  the  edge  before  being  de¬ 
formed  is  indicated  by  the  arrow. 

(b)  The  variation  of  the  real  area  of  contact  of  the  lead 
azide  crystal  shown  in  l  ig.  3a  when  it  is  compressed  by 
different  loads.  The  figure  on  the  left  of  each  frame  is  the 
load  (in  g)  on  the  crystal.  The  dark  band  represents  the 
real  contact,  whereas  the  bright  band  is  due  to  the  reflec¬ 
tion  of  the  light  front  the  plastically  deformed  surface 
when  it  is  only  slightly  separated  from  the  glass  plate. 

(c)  A  layer  of  particles  of  lead  azide  adhered  to  the 
glass  plate. 

increases,  both  the  width  and  the  length  of  the 
band  increase.  For  loads  above  24  g,  a  thin  bright 
line  appears  at  the  centre  of  the  dark  band.  At 
present  its  origin  is  not  clear.  The  crystal  was  then 
finally  unloaded,  and  during  the  unloading  the 
area  of  real  contact  decreases,  as  is  indicated  by 
the  change  in  the  contrast.  However,  as  the  crystal 
has  deformed  plastically,  the  decrease  is  slight 
initially  and  only  when  the  crystal  is  nearly  com¬ 
pletely  separated  from  the  glass  at  zero  load 
does  a  pronounced  decrease  in  the  area  of  real 
contact  occur.  As  already  stated,  the  flat  on  the 
crystal  forms  due  to  its  plastic  deformation;  it  is 
rather  smooth,  and  can  give  enhanced  illumination 
only  when  just  separated  from  the  glass  plate.  The 
behaviour  of  the  second  crystal  was  found  to  be 
quite  similar  to  that  of  the  first. 

It  should  be  pointed  out  here  that  if  a  crystal 
of  lead  azide  was  loaded  for  the  first  time  against 
the  glass  surface  even  by  as  low  a  load  as  I  g,  it 
left  on  the  glass  surface  a  thin  film  consisting  of 
very  small  particles  which  could  not  be  resolved 


(b) 

50  jim 


(a) 


/ 


(b)  0-1  mm 

Figure  4  (a)  A  crystal  of  LAT  type  C  before  being  com¬ 
pressed.  The  edge  to  lie  loaded  is  of  a  pyramidal  shape, 
with  the  angles  between  the  pairs  of  opposite  faces  being 
36°  and  81°.  (b)  After  the  compression  experiment  the 
tip  of  the  crystal  was  found  to  be  deformed  plastically 
(see  at  arrow),  (c)  The  variation  of  the  real  area  of  the 
contact  of  the  LAT  crystal  shown  in  fig.  4a  when  it  is 
pressed  against  the  glass  plate  at  different  loads.  The  load 
(in  g)  on  the  crystal  is  shown  on  the  left  of  each  frame. 
Dark  areas  represents  real  contact. 

even  under  a  high  powered  optical  microscope. 
An  examination  of  the  film  in  a  scanning  electron 
microscope  was  not  successful  either  because  of 
the  decomposition  of  the  particles  by  the  electron 
beam.  These  fine  particles  are  probably  of  lead 
azide  and  an  optical  micrograph  of  such  a  thin  film 
is  shown  in  Fig.  3c. 

Another  interesting  observation  was  that 
neither  during  the  approach  of  a  lead  azide  crystal 
towards  the  glass  plate  nor  during  the  crystal's 
removal  from  the  plate  after  having  been  loaded, 
were  the  surface  or  adhesive  forces  strong  enough 
to  be  registered  by  the  load  cell.  That  is,  these 
forces  were  less  than  0.25  g. 

3.2,  LAT 

We  examined  three  crystals  of  LAT  type  “C”  and 


«  « 

A 


(c)  100  jim 


only  one  of  LAT  type  “A".  The  compression 
behaviour  of  type  “C”  crystals  was  quite  similar 
and  results  from  a  crystal  are  shown  in  Fig.  4c. 
The  end  of  the  crystal  which  was  loaded  was  of 
pyramidal  shape  with  the  angles  between  the 
opposite  faces  being  36°  and  81°.  Figs.  4a  and  b 
show  this  end  before  and  after  deformation, 
respectively.  The  region  which  was  loaded  is  indi- 


(a)  0-1  mm 


cated  with  an  arrow  and  it  is  quite  clear  that 
plastic  deformation  of  the  tip  has  taken  place. 

It  will  be  seen  from  Fig.  4c  that  as  the  normal 
load  on  the  crystal  is  increased  from  0  to  Jig, 
the  area  of  real  contact  also  increases;  in  this  case, 
because  of  the  geometry  of  the  loaded  zone,  the 
area  of  real  contact  is  approximately  circular. 
When  unloading  occurs  and  the  load  is  reduced  to 
2g,  the  area  of  real  contact  only  reduces  slightly, 
as  in  the  case  of  lead  azide.  When  the  load,  is 
further  reduced  to  0.25  g,  some  separation  of  the 
surfaces  occurs,  as  is  evident  from  the  very  bright 
spot  in  the  middle  of  the  area  of  contact.  When 
the  surfaces  are  completely  separated  (i.e.  last 
frame),  there  are  still  some  dark  areas  left;  these 
are  due  to  an  adhered  layer  of  LAT  on  the  glass 
surface. 

Figs.  5a  and  b  show  a  crystal  of  LAT  (type  A) 
before  being  compressed  and  the  variation  of  real 
area  of  contact  with  load,  respectively.  As  will  be 
seen  in  Fig.  5a,  the  corner  of  the  crystal  which 
was  loaded  was  not  very  sharp  with  the  result  that 
the  contact  occurs  at  several  isolated  areas  (see 
Fig.  5b).  Moreover,  as  for  the  above  two  materials, 
when  the  load  was  increased  from  6.5  g  to  3 1  g, 
the  area  of  real  contact  also  increased.  On  com- 
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Figure  5  (a)  A  crystal  of  LAT  type  A  before  being  com¬ 
pressed;  the  loading  point  is  marked  with  the  arrow, 
(b)  The  variation  of  the  area  of  real  contact  of  the  LAT 
type  A  shown  in  l:ig.  5a  when  it  is  compressed  against 
the  glass  plate  at  different  loads.  The  load  (in  g)  on  the 
crystal  is  shown  on  the  left  of  each  frame.  Dark  areas 
represent  real  contact. 

plete  unloading,  a  very  thin  layer  of  adhered  par¬ 
ticles  (smaller  than  1  pm)  was  left  on  the  glass 
plate. 

3.3.  /3-HMX 

Two  single  crystals  of  this  material  were  investi¬ 
gated,  but  photographs  of  the  real  area  of  con¬ 
tact  of  only  one  of  them  will  be  shown  here. 
Fig.  6a  shows  the  face  of  the  crystal  that  was 
loaded  (this  crystal  did  not  have  any  pronounced 
corners  or  sharp  edges)  and  the  variation  of  the 
area  of  real  contact  with  load  is  shown  in  Fig.  6b; 
loads  in  grammes  are  given  at  the  bottom  left- 
hand  corner  of  each  frame.  In  this  case,  the  real 
area  of  contact  is  not  of  a  regular  shape.  Moreover, 
for  loads  greater  than  21.5  g  some  fringes  appear 
outside  the  area  of  real  contact.  Again,  for  these 
crystals  the  size  of  the  area  of  real  contact 
increased  with  load. 

The  plastic  deformation  of  a  sharp  corner  of 
the  second  crystal  is  clearly  illustrated  in  Fig.  7. 
Fig.  7a  shows  the  corner  of  the  crystal  which  was 
loaded  along  the  n-axis.  Fig.  7b  shows  the  state  of 
the  corner  after  the  crystal  had  been  loaded  to 
1 1  g  and  then  unloaded.  Note  the  permanent 
deformation  of  the  corner,  but  no  cracking  is 
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apparent  yet.  When  the  crystal  was  further  loaded 
to  23  g  and  then  unloaded,  the  volume  in  which 
permanent  deformation  occurred  also  increased 
(see  Fig.  7c).  This  frame  also  indicates  (hat  some 
fracturing  of  the  tip  may  have  occurred.  However, 
it  cannot  be  said  with  any  certainty  whether  the 
adhesion  of  the  crystal  tip  to  the  glass  enhanced 
the  fracture. 

It  may  he  added  here  that  a  film  of  adhered 
fine  particles  of  0-IIMX  crystals  was  left  on  the 
glass  surface  even  after  loading  to  only  3  g.  Whether 
this  film  formed  at  lower  loads  was  not  examined. 

3.4.  RDX 

Single  crystals  of  this  material  were  approximately 
spherical  and  of  diameter  less  than  I  mm.  Two 
crystals  were  investigated,  but  photographs  of  the 
area  of  teal  contact  of  only  one  of  them  will  he 
shown  here.  In  Fig.  8a  the  profile  of  a  crystal  is 
shown  with  an  arrow  marking  the  region  which 
was  loaded.  Fig.  8c  shows  the  variation  of  the  size 
of  the  area  of  real  contact  with  increasing  load; 
the  load  in  grammes  is  shown  a*  the  bottom  left- 


Figure  7  A  rather  sharp-edged  crystal  of  0-ltMX.  It  was 
first  compressed  to  a  load  of  1 1  g  (Fig.  7b)  and  then  to 
23  g  (Fig.  7c).  In  both  cases,  plastic  deformation  of  the 
crystal  tip  occurred. 

hand  corner  of  each  photograph.  At  zero  load 
the  crystal  is  very  close  to  the  glass  surface,  but 
not  in  contact  with  it.  This  causes  a  reflection  of 
the  light  as  is  indicated  by  the  lighter  area  in  the 
frame.  As  the  load  is  increased  to  1 .5  g  real  con¬ 
tact  occurs  at  the  point  indicated  by  the  arrow 
and  the  intensity  of  the  reflected  light  decreases 
giving  a  darker  zone.  Because  of  the  rough  nature 
of  the  crystal  surface  the  contact  occurs  at  several 
isolated  zones  as  is  clearly  indicated  by  frames 
corresponding  to  higher  loads.  This  sequence  gives 
a  clear  example  of  the  asperity  contact  described 
by  Bowden  and  Tabor  |13|.  The  maximum  load 
applied  to  the  crystal  was  50 g  (the  corresponding 
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Figure  8  (a)  Part  of  an  RI)X  crystal  which  did  not  exhibit 
any  sharp  edges.  It  was  loaded  at  the  region  marked  with 
the  arrow,  lb)  After  loading  to  41  g,  a  thin  sliver  of  the 
crystal  appears  to  be  on  the  point  of  detaching  when  the 
crystal  was  separated  from  the  glass  plate. 
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f  igure  9  The  variation  nf  (tv  area  of  real  contact  of  a  wedge-like  crystal  of  KCI  when  it  was  compressed  along  1 1  1  0 ) 
against  the  glass  surface  at  dillerent  loads.  The  load  (in  g)  is  shown  at  the  bottom  left-hand  corner  of  each  frame.  Dark 
areas  represent  the  real  contact.  Note  also  the  interference  fringes;  these  correspond  to  regions  outside  (he  real  contact. 


area  of  real  contact  is  not  shown  here)  and  then 
it  was  unloaded  and  again  a  very  thin  layer  of 
extremely  small  particles  was  left  adhering  to  the 
glass  plate.  A  photograph  of  the  crystal  after  the 
deformation  is  shown  in  Tig.  8b.  It  appears  that  a 
thin  sliver  of  the  crystal  was  in  the  process  of 
being  torn  off  at  the  contact  zone. 

3.5.  KCI 

F:ig.  9  shows  the  variation  of  the  real  area  of  con¬ 
tact  for  a  single  crystal  of  KCI  loaded  on  a  fresh 


edge  formed  by  cleaving  along  {100}  planes;  the 
direction  of  compression  was  approximately  along 
1 1  1  0|.  (Note  that  it  is  because  of  the  lighting  con¬ 
ditions  that  in  (he  figure  only  one  of  (lie  {1  00} 
faces  forming  the  edge  is  seen.)  In  the  first  frame, 
the  load  is  2.5  g  and  real  contact  occurs  at  isolated 
points  marked  by  the  arrows;  interference  fringes 
are  also  clearly  visible  outside  the  real  contact 
zone.  As  the  load  is  increased,  the  number  of 
zones  in  real  contact  and  their  areas  increase,  but 
no  cracking  of  the  crystal  occurs.  The  load  was 
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increased  to  25.5 g  and  then  the  crystal  was  un¬ 
loaded.  Note  that  until  the  load  has  dropped  to 
zero  there  is  no  significant  reduction  in  the  area 
of  real  contact  (see  the  frame  corresponding  to 
7.5  g).  This  strongly  suggests  that  real  area  of 
contact  is  formed  by  plastically  deformed  area  of 
the  material.  Naturally,  as  the  load  is  decreased, 
the  area  of  the  plastically  deformed  flat  will  not 
reduce.  When  the  load  is  reduced  to  zero  and  the 
surfaces  are  just  separated,  the  plastically  flowed 
material  gives  enhanced  illumination.  In  the  last 
frame  (i.e.  zero  load)  a  very  small  zone  is  still  in 
contact  (dark  zone  marked  with  arrow);  this  dark 
zone  also  changed  to  a  bright  one  on  further 
separation.  For  this  crystal  no  layer  of  transferred 
material  was  left  on  the  glass  plate. 

4.  Data  reduction 

The  average  contact  stress  was  obtained  by  dividing 
the  applied  load  by  the  corresponding  measured 
area  of  real  contact.  As  in  all  cases,  even  at  the 
smallest  loads  studied,  plastic  deformation  of  the 
material  had  occurred  over  the  entire  contact 
zone,  the  average  contact  stress  gives  us  the  plastic 
deformation  stress  of  the  material.  The  variation 
of  the  deformation  stress  with  load  of  the  differ¬ 
ent  crystals  are  shown  in  Figs.  10  to  14.  Except 
for  KC1,  all  crystals  studied  show  that  the  defor¬ 
mation  stress  increases  with  load,  reaching  a 
plateau  at  higher  loads.  For  crystals  of  different 
materials  the  average  deformation  stress  at  high 
loads  was  greater  by  a  factor  in  the  range  1.5  to 
5.0  than  that  at  low  loads.  Notice,  however,  that 


Figure  10  The  variation  of  the  deformation  stress  of 
o-lead  azide  crystals  with  load.  •,  crystal  No.  I ;  *,  crystal 
No.  2. 


Figure  1 1  The  variation  of  the  deformation  stres*  of  LAT 
type  C  and  type  A  single  crystals  with  load.  •,  a,  c,  LAT 
type  C  crystal  Nos.  1,  2  and  3  respectively;  a,  LAT  type  A 
single  crystal. 

deformation  stress  against  load  graphs  show  irre- 
producible  behaviour  even  for  crystals  of  the  same 
material.  Several  factors  may  contribute  towards 
this  irreproducibility:  (i)  a  change  of  the  loading 
axis  of  the  crystal;  (ii)  a  change  in  the  geometry 
of  the  loaded  region;  and  (iii)  changes  in  the 
physical  properties  from  crystal  to  crystal  due  to 
the  impurity  content  variations. 

For  the  case  of  a-lead  azide  results  from  exper¬ 
iments  on  two  crystals  are  shown  (Fig.  10).  For 
one  crystal  the  deformation  stress  increases  from 
13xl07Pa  to  26  x  107Pa  as  the  load  increases 
from  about  10-’  N  to  4  x  10"'  N,  whereas  for  the 
second  crystal  the  deformation  stress  increases 
from  5xl07Pa  to  10xl07Pa  with  the  increase 
of  load  from  4  x  10_,N  to  2  x  10'1  N.  The  results 
for  three  crystals  of  LAT  type  C  and  one  of  LAT 
type  A  are  shown  in  Fig.  11.  Crystal  No.  1  shows 
an  increase  of  deformation  stress  of  over  four 
times  (i.e.,  from  ~5xl07Pa  to  22  x  107Pa)  as 
the  load  goes  up  from  2  x  10_JN  to  40  x  10_1N. 
However,  in  the  other  two  crystals  the  variation  of 
the  deformation  stress  is  not  so  marked;  crystal 
No.  2  shows  an  increase  of  only  25%  in  the  defor¬ 
mation  stress,  and  from  the  small  number  of  data 
points  (four  only)  covering  a  relatively  small  load 
range  (4  to  10  x  10*JN)  for  crystal  No.  3,  very 
little  variation  seems  to  occur. 

In  d-HMX  crystals  (Fig.  12)  the  variation  of  the 
deformation  stress  with  load  is  from  25%  (crystal 
No.  1 )  to  a  factor  of  2  (crystal  No.  2),  though  the 
load  range  of  crystal  No.  2  is  considerably  larger. 
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Figure  12  The  variation  of  the  deforma¬ 
tion  stress  of  0-HMX  single  crystals  with 
load.  crystal  No.  1 ;  •,  crystal  No.  2. 


One  significant  difference  is  the  value  of  the  stress 
at  low  loads;  the  deformation  stress  of  one  crystal 
is  3  to  4  times  that  of  the  other.  In  this  case  it 
does  appear  that  the  geometry  of  the  contact  is 
playing  an  important  part;  the  lower  stress  values 
are  for  the  slightly  rounded  crystal. 

RDX  crystals  (Fig.  13)  also  show  an  increase 
in  the  value  of  the  deformation  stress  by  a  factor 
of  2,  as  the  load  is  increased  from  2  x  10~J  N  to 
40  x  l  O'*  N.  However,  as  for  the  other  materials, 
there  is  a  difference  in  the  stress  values  of  one 
crystal  from  that  of  the  other. 

For  the  KC1  crystal  (Fig.  14)  the  deformation 
stress  remains  constant  at  2.7  x  107Pa  as  the  load 
increases  from  4  x  10'J  N  to  25  x  10"J  N. 

It  may  be  pointed  out  that  for  all  the  crystals 
studied,  the  size  of  the  plastically  deformed  zone 
at  the  contact  was  not  more  than  l/50th  of  the 
size  of  the  crystal.  In  other  words,  the  plastic 
deformation  was  very  localized. 


Figure  13  The  variation  of  the  deformation  stress  of  RDX 
single  crystals  with  load,  a,  crystal  No.  1 ;  a,  crystal  No.  2. 


5.  Discussion 

From  the  analysis  of  the  data  given  in  Section  4, 
several  important  observations  can  be  made  regard¬ 
ing  the  flow  stresses  of  small  (mm-size)  single 
crystals  of  a  number  of  brittle  materials: 

1 .  plastic  flow  occurs  even  for  loads  as  small  as 
0.02  N; 

2.  generally,  the  flow  stress  increases  with  the 
applied  load  reaching  a  plateau  at  loads  of  10_l  to 
2x10'*  N; 

3.  the  maximum  deformation  stress  of  a  crystal 
is  considerably  smaller  than  its  Vickers  diamond 
hardness  value  (see  Table  II); 

4.  except  for  KC1,  all  crystals  studied  leave  a 
film  (composed  of  extremely  small  particles)  on 
the  glass  surface  against  which  they  are  com¬ 
pressed,  even  for  loads  as  small  as  10'J  N. 

We  shall  discuss  these  observations  separately 

5.1.  Plastic  flow  at  very  low  loads 

The  crystals  used  were  as-grown  and  it  is  very  likely 


Figure  14  The  variation  of  the  deformation  stress  of  a  KC1 
single  crystal  with  load. 
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that  their  mechanical  behaviour  would  be  that  of 
non-work-hardened  crystals.  In  fact,  using  the 
technique  of  indentation  with  a  hard  sphere,  it  has 
been  found  that  the  crystals  of  0-HMX  and  RDX 
show  a  considerable  amount  of  work  hardening. 
Now  since  some  of  the  crystals  possess  “point¬ 
like”  corners  and  some  have  “line-like”  edges,  the 
elastic  limit  will  be  exceeded  at  extremely  low 
loads.  This  has  been  confirmed  by  the  experimental 
observations  reported  in  Section  3. 

5.2.  The  increase  of  the  deformation  stress 
with  load 

The  increase  of  the  deformation  stress  with  load 
is  a  little  difficult  to  explain  fully  for  conical, 
pyramidal  and  wedge-shaped  crystals.  It  may  be 
possible  that  a  very  limited  number  of  slip  systems 
available  in  these  crystals  is  a  contributory  factor. 
It  should  be  pointed  out  that  in  the  case  of  cones 
and  wedges  of  an  isotropic  material  possessing 
more  than  five  independent  slip  systems,  such  as  in 
copper,  the  deformation  stress  is  independent  of 
the  load  for  the  fully  work-hardened  state  and 
only  increases  by  about  6%  when  the  copper  is 
annealed  (14],  Therefore,  for  the  brittle  crystals 
studied,  tt  is  reasonable  to  assume  that  the  work 
hardening  effects  will  show  some  increase  in  the 
deformation  stress  with  the  increase  of  load. 
However,  more  experimental  work  is  needed  on 
other  brittle  crystals  before  almost  a  100%  increase 
in  the  deformation  stress  (see  Figs.  10  to  13)  with 
load  can  be  understood. 

For  crystals  whose  loading  edges  were  curvedfe.g. 
RDX),  the  deformation  stress  against  load  behav¬ 
iour  is  somewhat  similar  to  that  of  an  annealed 
metallic  sphere;  in  both  cases  the  stress  initially 
increases  with  load,  reaching  a  peak  at  higher  loads. 
On  comparison  with  the  results  from  metallic 
spheres  ( 1 5 J,  it  appears  that  the  initial  yield  in  the 


crystal  occurs  at  a  value  only  just  a  few  per  cent 
higher  than  its  uniaxial  yield  stress.  The  increase 
of  the  deformation  stress  with  load  can  then  be 
explained  by  the  work  hardening  of  the  crystal. 

5.3.  The  maximum  deformation  stress  of 
a  crystal  against  its  Vickers  diamond 
hardness 

The  maximum  deformation  stress  and  Vickers 
diamond  hardness  of  the  various  crystals  are  com¬ 
pared  in  Table  II.  It  is  noted  that  in  all  cases  the 
latter  is  considerably  higher.  Two  factors  respon¬ 
sible  for  the  difference  are: 

(i)  the  constraint  imposed  on  a  solid  when  a 
Vickers  diamond  pyramid  is  pushed  into  the  solid 
is  greater  than  when  a  corner  or  an  edge  of  the 
solid  is  pressed  against  a  hard  flat  surface;  and 

(ii)  the  average  plastic  strain  induced  by  the 
Vickers  diamond  pyramid  indentation  is  8%, 
whereas  the  strain  induced  due  to  the  plastic 
deformation  of  a  sharp-edged  or  pointed  particle 
is  dependent  upon  the  included  angle  of  the  loaded 
zone,  as  has  been  shown  for  the  case  of  annealed 
copper  cones  1 14). 

One  important  implication  of  these  findings  is 
that  to  predict  the  real  area  of  contact  of  a  brittle 
particle  loaded  against  a  hard  flat  surface  using  its 
Vickers  diamond  hardness  value,  as  has  been  done 
by  the  previous  investigators,  is  not  generally 
correct;  this  approach  gives  very  much  an  under¬ 
estimate. 

5.4.  The  layer  left  on  the  glass  plate 
Except  for  KC1,  all  crystals  left  a  layer  of  micro¬ 
scopic  particles  on  the  glass  plate  even  when 
compressed  by  as  low  a  load  as  1  g.  There  appear 
to  be  two  possible  explanations: 

(i)  the  crystals  break  up  due  to  the  compression; 
and 


TABLE  II  Comparison  of  maximum  deformation  pressure  of  a  particle  and  its  Vickers  diamond  hardness 


Material 

Vickers  diamond 
hardness*  //v  (kg mm'1) 

Maximum  deformation 
pressure  of  particle 

Pm  (kgmm1)^ 

fm/tfv 

119  .7  *  8.1 

27.5 

0.23 

(PbN,) 

RDX 

24  1  ?  0.8 

15 

0.62 

0-H  M  X 

41.3  *  10 

26 

0.63 

LAT  (type  C) 

109  4  t  3,1 

29 

0.27 

LAT  (type  A) 

34  5  t  0  9 

12 

0.35 

KC1 

16  !  1 

2.8 

0.18 

“Load  on  inrientor  was  1 5  g 
+ Error  =»  5'.(. 
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TABLE  111  Predicted  minimum  loads  for  compression 
cracking 


Materials 

Predicted  minimum* 

(all  single  crystals) 

load  for  cracking  (N) 

a-lead  azide 

6.6  X  10“ 

RDX 

1.2  x  10“ 

KC1 

1.9  X  10“ 

’Calculated  using  //y  values. 

(ii)  the  break  up  does  not  occur  during  the 
compression,  but  the  adhesion  between  the  crystal 
and  the  glass  surface  is  so  high  that  when  the 
crystal  is  separated  from  the  glass,  fracture  occurs 
within  the  crystal. 

According  to  Hagan  [16],  the  critical  load  Pa 
to  cause  cracking  of  a  brittle  particle  by  com¬ 
pression  is  given  by: 

Pa  =  m{KlcIHv)3Klc  (3) 

where  KIC  is  the  critical  stress  intensity  factor  of 
the  particle  and  tfv  its  Vickers  diamond  pyramid 
hardness.  For  plane  strain  conditions,  Ktc  = 
[(2yE)/(l  —  v2)]1'2  where  y  is  the  fracture  surface 
energy,  E  the  Young’s  modulus,  and  v  the  Poisson’s 
ratio.  The  critical  calculated  loads  to  cause  crack¬ 
ing  are  shown  in  Table  III. 

It  will  be  seen  from  Table  III  that  the  crystals 
of  PoN6,  RDX  and  KC1  should  crack  at  the  loads 
used  in  these  experiments.  It  was  pointed  out 
above  that  the  deformation  pressures  of  the  various 
crystals  are  considerably  smaller  than  their  Vickers 
hardness  values  (Table  II).  Therefore,  if  we  use  the 
maximum  deformation  pressures  in  Equation  3  in 
place  of  Hv,  the  predicted  loads  will  be  markedly 
higher,  though  both  lead  azide  and  RDX  will  still 
crack,  but  not  KC1,  for  the  loads  employed  here. 

It  must  be  pointed  out  that  in  situ  observations 
of  the  explosive  crystals  did  not  show  any  crack¬ 
ing,  suggesting  that  the  adhered  layer  on  the  glass 
probably  forms  during  the  unloading  (i.e.  separ¬ 
ation).  However,  irrespective  of  the  fact  whether 
cracking  occurs  during  loading  or  unloading,  the 
measurement  of  the  area  of  real  contact  will  not 
be  affected. 

The  fact  that  when  a  particle  of  an  explosive 
material  comes  into  contact  with  a  solid  surface 
(e.g.  a  container),  it  leaves  a  thin  layer  of  the 
explosive  on  the  solid,  means  that  the  situation 
can  be  hazardous.  Therefore,  it  is  important  to 
clean  the  surface  of  the  container  thoroughly  after 
it  has  been  in  contact  with  an  explosive,  especially 
if  it  is  of  a  primary  type. 


5.5.  Flow  stress  determination 
From  the  deformation  stress  values  of  a  brittle  par¬ 
ticle,  an  estimate  of  its  flow  stress  (i.e.  yield  stress 
at  the  imposed  strain)  can  be  made.  If  the  particle 
is  spherical,  then  the  deformation  stress  at  low 
loads  will  correspond  closely  to  the  uniaxial  yield 
stress  of  the  explosive  crystals  used  in  this  work. 
The  values  are: 

a-PbN6  <  5  x  107  Pa; 

LAT  type  C  <  5  x  107  Pa; 

LAT  type  A  <  6.5  x  107  Pa; 

(3-HMX  <  7.5  x  107  Pa; 

RDX  <  4  x  107  Pa. 

6.  Conclusions 

The  plastic  deformation  stress  of  as-grown  single 
crystals  of  a  number  of  explosive  materials  has 
been  successfully  determined  using  a  new  exper¬ 
imental  technique.  The  experiments  involved 
measuring  the  real  area  of  contact  between  the  test 
crystal  and  a  transparent  glass  plate  pressed 
together  with  a  given  normal  load.  The  deforma¬ 
tion  stress  was  found  to  increase  with  load,  reach¬ 
ing  a  plateau  at  higher  loads.  It  was  discussed 
whether  the  plastic  deformation  stress  of  a  particle 
approaches  its  uniaxial  yield  stress  at  low  normal 
loads;  estimates  of  the  limiting  values  of  the  uni¬ 
axial  yield  stress  of  the  explosive  crystals  have 
been  given. 
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